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ON COMPATIBLE DIFFERENTIAL EQUATIONS 
AND THE ORTHOGONAL PROPERTIES OF 
THEIR SOLUTIONS 


By W. 0. KERMACK (Edinburgh) AND W. H. McCREA (London) 
[Received 2 December 1932] 


1. Introduction. In some previous work} we have defined con- 
jugate partial differential equations for a function of two variables. 
In the present paper we begin by examining more fully the condition 
for their compatibility and that of their adjoint equations. We then 
extend the work to the case where, instead of partial derivatives with 
respect to one of the variables, we have differences with respect to 
that variable. This is the case applicable to discrete eigenvalues. It 
is well known that the eigensolutions of a linear differential equation 
involving a parameter are orthogonal to those of its adjoint equation. 
This parameter can be treated as the second variable in the present 
work, and a difference equation with respect to its eigenvalues can 
be written down for the eigenfunctions. What we now show is that 
the eigensolutions of this differential and difference equation, and 
those of the adjoint pair of equations, form normalized orthogonal 
sets. 

We do not examine fully the general conditions under which the 
latter theorem is true. Our purpose is rather to show how it can be 
used to predict the existence and form of orthogonality relations, 
the strict validity of which may require independent investigation. 
Treated in this way the work brings out a fresh aspect of the relation 
between the differential and difference equations satisfied by func- 
tions of hypergeometric type, and the relation of these equations to 
the orthogonal property. 

Although the present paper is more or less self-contained, we retain 
the notation of previous papers, since we hope elsewhere to develop 
the theory from the point of view of the algebra of operators. 


2. Conjugate differential equations. We denote by g, Q two 
independent variables, and by p, P the operators —é/éq, 0/eQ 

+ W. O. Kermack and W. H. McCrea, Proc. Edinburgh Math. Soc. (2) 2 
(1931), 205-19, 220-39; Proc. Roy. Soc. Edinburgh, 41 (1931), 176-89. These 
will be referred to as Papers I, II, III. Conjugate equations are defined in 
Paper II. 
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respectively. A pair of partial differential equations for the function 
x(q, Q) of the form 


[Q(¢.p)—@|x = 0 (2.01) 
[P(q.p)—P]x = 9, (2.02) 


where Q(q,p), P(g, p) are given functions of g,p, will be called con- 
jugate equations in normal form. We have 


THEOREM I. A necessary and sufficient condition that the equations 

(2.01), (2.02) should be compatible is that 
Q(g. P)P.p)—Pg. PAG P) = 1 (2.03) 
identically, where qp—pq = 1. 

This was given as Theorem XIII of Paper II. In order, however, 
to make explicit the degree of arbitrariness of the common solution 
x(q, Q), we shall give an independent proof of the sufficiency of the 
condition (2.03). This amounts, as a matter of fact, to finding of how 
many independent interpretations the function w(q,p) of Paper II 
is susceptible, when treated as an operator. We hope later to con- 
sider this problem directly. We now show that 

THEOREM II. Jf the condition (2.03) is satisfied, and if Q(q,p) is 
an integral function of p of degree n, then there are in general n, and 
only n, common solutions x(q, p) of (2.01), (2.02), linearly independent 
in both q and Q. 

We shall write out the proof for n = 3, but the method is general. 

Equation (2.01) by itself is a linear ordinary differential equation 
for x as a function of g, and will possess three solutions linearly 
independent in g. These we denote by y,(¢, Q), x2(¢, 2), x3(q; Y); 80 
that the general solution of (2.01) is 


X = % X11 X2+4%s Xz; (2.04) 
where a, @,, a; are arbitrary functions of Q. 
Since [Q(¢.P)—@lx1 = 9, (2.05) 
then [P(¢. p)Q(q,P)—P(G. P)@]x1 = 0 
or, using (2.03), 
(Q(q.p)P(q, p)—@P(G.p)—1]x1 = 0. (2.06) 


Similarly we have again, from (2.05), 
[PQ(¢,p)—P@]x = 0 
or, since PQ—QP = 1, 
[Q(q. p)P—QP—1]x, = 0. (2.07) 
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Hence, subtracting (2.06), (2.07), we have 


[Q(q.p)—QIP@. p)—P]x = 0. (2.08) 
Therefore [P(q, p)— P }x, is itself a solution of (2.01), and so must be 
of the form (2.04). This we may express, together with the corre- 
sponding results for x5, x3, by saying that functions a,,, @,»,..., G33 of 
Q exist, such that 


[P(q. p)—P]xa = 941 X1 +412 X2+s Xs° (2.09) 
[P(q, p)—P]x2 = 41 X1 +422 X2+4123 Xa» (2.10) 
[P(q, p)—P]xs = 431 X1 +452 X2+4538 X3- (2.11) 


Multiplying these equations by functions b,(Q), b.(Q), 63(Q) re- 
spectively, and adding, we obtain 
[P(q,p)—P]x = (ay, 6; —; +4; 62+45; b3)x4 + 


+ (Ay2b1+-G 92 bg—b3 +39 b3)x2+ 


+ (443 by +493 6, +4133 63 —b3) x3, (2.12) 
where b; = db,/dQ, etc., and where 
X = by x1 +52 x2+5s xz. (2.13) 


Since x;, Xa, X3 are linearly independent, it follows from (2.12) that 
x is a solution of (2.02), that is, 

[P(q.p)—P]x = 9, 
if and only if the coefficients of x,, x2, x3 on the right-hand side of 
(2.12) vanish separately, giving 


a1, 6, —b| +g, by +a; bs = 0, (2.14) 
yo 6, +9 b,—b3; +55, = 0, (2.15) 
43 6, +Ao3 b.+-A336,—6, = 0. (2.16) 


These equations determine the functions 6,, b,, b;. The elimination 
of any two of them yields in general a linear differential equation of 
order three in the remaining one. To each of its three solutions there 
will then correspond one set of functions b,, b., bs. Hence, in general, 
there are three, and only three, functions of the type x, and the 
theorem is proved. Any linear combination of these functions with 
constant coefficients will be a common solution of (2.01) and (2.02). 


3. Adjoint equations. We define the adjoint f*(q,p) of a func- 
tion f(q, p) of g,p as the function obtained by reversing the order of 
the factors in each term of f(g,p), and at the same time changing 
the sign of p. Thus, for example, if f(qg,p) = q’p?+qp+1, then 

G2 
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f*(q, p) = p?g?—pq+1. In particular, the adjoint of a product of two 
functions f(q, p), g(¢, p) is given by 

(f9)* = 9*f*. 

The adjoint of a given differential equation is then defined as the 
equation obtained from it by replacing all the functions of operators 
occurring by the adjoint functions. 

Now it has been shown by N. H. McCoy7 that the adjoint of an 
identity is itself an identity. Applying this to the identity (2.03), 
we Rave P*(q, p)Q*(q, p)—Q*(q, p)P*(q, p) = 1. (3.1) 
But from Theorem I this is precisely the necessary and sufficient 
condition that the pair of differential equations 

[Q*(q,p)—Q]x* = 0 (3.2) 

[P*(q, p)+P]x* = 0 (3.3) 
for the function y*(q,Q) should be compatible. Now these are the 
equations adjoint to (2.01), (2.02). Thus we have shown 

THEOREM III. The adjoint equations of a pair of compatible dif- 
ferential equations in normal form are themselves compatible and in 
normal form. 

Clearly the number of independent solutions of (3.2), (3.3) is given 
by Theorem II, and is the same as the number of solutions of (2.01), 
(2.02). For if Q(q, p) is of degree n in p, then Q*(q, p) is also of degree 
nin p. 

4. Second type of compatible equations. We shall be led by 
later applications to consider a pair of equations for the function 
x(q, Q) of the form [Qq.p)—@lx = 0 (4.1) 

[S(q, p)—e? ]x = 0, (4.2) 
where Q(q, p), S(q, p) are given functions of g,p, and where we write 
formally for any function of Q, f(Q+1) = e’f(Q). Hence (4.2) is 
actually a difference equation for y. 

The following results can now be proved by methods exactly 
similar to those already used for the corresponding theorems in §§ 3, 4. 

THEOREM IV. A necessary and sufficient condition that equations 
(4.1), (4.2) should be compatible is 

Q¢. PS. P)—-S(Z. P)QY. Pp) = Sq. P) (4.3) 
identically, where qp—pq = 1. 
+ N.H. McCoy, Trans. American Math. Soc. 31 (1929), 793-806, Theorem IT. 
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THEOREM V. If the identity (4.3) is satisfied, and if Q(q,p) is an 

integral function of p of degree n, then there are n, and only n, common 
solutions of (4.1), (4.2), linearly independent in both q and Q. 

The proof of Theorem V follows the steps of Theorem I, but leads 
to a difference equation of order n for each of the functions b,, bg,.... 
The n linearly independent solutions of this difference equation can 
be multiplied by arbitrary functions of Q of unit period. But since 
we shall in all applications confine attention to values of Q at unit 
interval, these periodic functions have no effect, and can be ignored 
from the outset. 

THEOREM VI. The identity (4.3) gives also a necessary and sufficient 
condition that the pair of equations adjoint to (4.1), (4.2) should be 
compatible. 

The adjoint equations are 

[Q*(q, p)—Q]x* = 0 (4.4) 
[S*(q, p)—e-* ]x* = 0. (4.5) 
If they are compatible they possess » solutions linearly independent 
in both g and Q. 
5. Orthogonal properties. We shall now show 


THEOREM VII. Jf a pair of compatible difference-differential equa- 
tions are given in the form 


[Q(q,p)—@]x = 9, (5.01) 
[S(q,p)—e*]x = 0, (5.02) 
and their adjoint equations in the form 
[Q*(q, p)—Q]x* = 9, (5.03) 
[S*(q, p)—e-P ]x* = 0, (5.04) 


and if a pair of solutions x(q, Q), x*(q, Q) can be made to satisfy suitable 
boundary conditions at end-points q = a, q = b (say), by restricting Q 
to assume integral eigenvalues, then the functions x, x* are normalized 
orthogonal functions with respect to the interval (a,b), i.e 
b 
| x(q, Q)x*(q, @’) SQ" (5.05) 
a 
where 899: = 0, 9 4 O'; 899 = 1, and where k is a purely numerical 
constant independent of Q. 
Suppose Q(qg,p) is expanded as a sum of terms of the form q‘p/ 
where, since we still assume Q(g, p) to be an integral function of p, 
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j must be a positive integer. Consider, then, the integral 
b 

| {a'pix}x* dg 


a 








= [—p!y.q'x*—pl*x. (—p)q'x*—...—x- (=p gx + 
b 
+ | x{(—p)a'x*} dg (5.06) 
b a 
= [...2+ | xt(qip!)*x*} dq (5.07) 


by elementary integration by parts. Treating each term of Q(q, p) 
in this way and adding the results, we have, using (5.01) and (5.03), 
b b 
Q | x(q, Q)x*(q, Q') da = | (Qa. P)x(q. Q}x*(q. Q) 


a a 


b 
= | x(q, Q4Q*(@.P)x*(q Q')} dg 


a 
b 


= Q' | x(q. @x*(a. @') da, (5.08) 
a 
together with an integrated part, given by the sum of terms corre- 
sponding to those in square brackets in (5.06). 

The condition that this integrated part should vanish for the end- 
points a, b, for all permitted values of Q, Q’ provides the eigenvalue 
problem. This may or may not have a solution. _ In a certain class 
of cases it is solved by taking integral eigenvalues for Q; and we 
confine ourselves for the moment to such cases. 

Hence we have from (5.08) for the corresponding eigenfunctions 


(Q—Q’)% = 0, where 
£(Q, Q’) = | x(q, Q)x*(¢. Q’) dq. (5.09) 


Using equations (5.02), (5.04), and proceeding as in the derivation 
of (5.08), we have similarly 
b b 
eP | x(q, Q)x*(q, Q') dq = e-® | x(q, Q)x*(q,@') dq, (5.10) 
together with an integrated part, where P’ is conjugate to Q’. Now 
we shall assume that we are dealing with cases where the same 
boundary conditions that cause the integrated part of (5.08) to vanish 
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cause also the integrated part of (5.10) to vanish. The examples we 
shall give serve to show that this holds in a large class of cases, but 
we shall not here investigate the general conditions of its validity. 
Hence from (5.10) 

(Pe) = 0, or HQ4+1,Q)—W(Q,Q’'—1) = 0. (6.11) 


Now (5.09) expresses the usual result for eigenfunctions, namely, 


¥(Q, Q') = K(Q)8oq°, (5.12) 
where K(Q) is some function of Q (provided the integral has no 
singularity at @ = Q’). 

Writing Q’ = Q+1 in (5.11) and using (5.12), we now have 
K(Q+1) = K(Q) = k, say, (5.13) 
where k is a constant independent of Q. This is the required result. 
The orthogonal property of the eigensolutions of a single linear 
differential equation and its adjoint equation are, of course, well 
known. The extension here introduced is in showing that the eigen- 
solutions of a compatible pair of equations, and the adjoint pair, are 
orthogonal and normalized. 


6. More general form of the equations. A more general form 
covering all the examples we shall consider is given by a pair of 


equations [f(a.p)—F(Q,e”)]x = 0, (6.01) 
[9(q, p)—G(Q, e”) |x = 0, (6.02) 
where f,g, F,@ are given functions, and where the operators have 


no common factors. 

We shall not investigate a sufficient condition for the compatibility 
of a pair of equations of this type. But we shall suppose that a 
common solution x(q, Q) is known to exist. Then if this x(q, Q) were 
given, we could by the methods described in our Paper I derive a pair 
of conjugate equations formally satisfied by x, and by a further 
application of the same formal methods we can transform these into 


the type [Q@7)—-Qlx = 9,  [S@p)—e”]x= 0. (6.08) 
Since, then, y is a solution of (6.01), (6.02), we have by substituting 
the values of Qy, e?x given by (6.03), the relations 
f(a. P)—F{Q@ P). Sq, p)} = 0 (6.04) 
9(7. P)-G{Qq, p), Sq. p)} = 0 (6.05) 
satisfied identically, where gp—pq = 1, and where F, G@ are the 
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transposed functions F, G obtained by reversing the order of the 
factors in every term of F, G. 
By McCoy’s theorem already quoted we have now from (6.04), 
(6.05) the adjoint identities 
f*(¢,p)—F*{Q*(q, p), S*(q, p)} = 0, (6.06) 
9*(q, p) -G*{Q*(q, p), 8*(q, p)} = 0. (6.07) 
Therefore any solution of the equation 
[Q*(q, p)—Q]x* = 0, [S*(q, p)—e-']x* = (6.08) 
is also a solution of the equations} 
[f*(q,p)—F*(Q, e?)]x* = 0, (6.09) 
[9*(q, p)—G*(Q, e?)]x* = 0. (6.10) 
But since equations (6.03) are by their derivation compatible and 
have the solution y, it follows by Theorem VI that equations (6.08) 
are compatible and have a common solution y*. Hence we have 
shown that the equations (6.09), (6.10) adjoint to equations (6.01), 
(6.02) are compatible, and have the common solution x*. Further, 
provided suitable boundary conditions can be satisfied by taking 
integral eigenvalues for Q, then by Theorem VII applied to (6.03), 
(6.08), the functions y, y* are orthogonal and normalized. So we have 
THEOREM VIII. Jf a pair of distinct difference-differential equations 
(f(a. p)—F(Q,e")|x = 9, = [9(¢,. p) —G(Q, e”) ]x = 0 


possess a common solution x(q, Q), then the adjoint equations 
[f*(q.p)—F*(Q,e")]x* = 0, — [g*(¢, p) —G@*(Q, e”)]x* = 
possess a common solution x*(q,Q), and, provided suitable boundary 


conditions at end-points q =a, q = b (say) are satisfied by taking 
integral eigenvalues for Q, we have further 

b 

[ x(a, Q)x*(a, @') dq = kB gq. 

a 


where k is a constant independent of Q. 


7. Non-integral eigenvalues. If the boundary-value problem 
leads to non-integral eigenvalues for Q, we can still apply the theory, 
provided the eigenvalues Q are expressible as functions of an integral 
variable N. For we can then change the variable from Q to N and 


+ Note that it is correct to write F*(Q,e”), and not F*(Q,e-”). For 
example, if F(Q, e”) = Qe”, then F*(Q, e”) = e~”Q gives the adjoint function. 
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equations of the form (6.01), (6.02) will transform into others of the 
same type. These will have integral eigenvalues and Theorem VIII 
applies. 


8. Continuous eigenvalues. A similar theory can be written 
down formally for the case of continuous eigenvalues, with Dirac’s 
5-function replacing the Kronecker 5-symbol. But the satisfaction 
of the boundary conditions appears to be more difficult, and we have 
not succeeded in constructing examples for which it holds. 


9. Applications and examples. As an example of a pair of con- 

jugate equations in normal form we may take the pair 
(P?+Q)x=0,  (dgp1—P)x = 0. 
Here we have Q(q, p) = —p*, P(q.p) = 4¢p-, so that 
Q(q. p)P(Q. P)— Pd. P)QG. P) = —3P*ap*+h9p Pp? = 1, 

and therefore the identity (2.03) is satisfied. 

We find by direct solution that there are two, and only two, 
linearly independent solutions 

x = singvQ, x = cosqvQ, 

illustrating the result of Theorem IT. 

The adjoint equations are 

(pP?+Q)x*=0, (3p-"g+P)x* = 0. 

We find again by direct solution that these possess just two 

independent solutions 


- 1 
= —singvQ, * = —.cosqvg@. 
. ae Q ir mek Q 


This verifies Theorem III. 


10. Legendre polynomials P,(z). The usual differential and 
difference equations for P,(z) are 


[a2 25 +n(n+1)|PyC) = 0 (10.1) 
dz* Z 

(n+I)P, 4(2)—(2n+ YjzP,(z)+-nP,_(2) = 0. (10.2) 

To adhere to our customary notation we shall write P,(z) as 

x(q; Q), where q, Q replace z, n respectively. Then (10.1), (10.2) be- 

77 [pA—9)p+Q(Q+1)]x = 0 (10.3) 


ee _ 
E 2041° , e : (10.4) 
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These equations are by their derivation compatible, and we can 
apply Theorem VIII. 

The adjoint equations are 


[p(l—q?)p+ Q(Q+ 1)]x* = (10.5) 
san» @ i, 
la—e F 3041 —el oa* = 0. (10.6) 


Since the first equation is self-adjoint and does not involve the 
operator P, we try a solution 
x*(q, Q) = B(Q)x(@, @), (10.7) 
where A(Q) is a function of Q only. We find that this gives a solution 
of (10.6) if 8 satisfies the equation 
B(Q+1)__ BQ) _ B(Q—1) 
20+3 2Q+1 2Q-1° 
That is, we may take B(Q) = 2Q+1. 
Thus the compatibility of the adjoint equations is verified. Theorem 
VIII will further give the orthogonal property 
b 
| x(q, Q)x*(q, 2’) dq = k8gq,, (10.8) 
provided suitable end-points a, b can be found. But it is clear from 
the form of the differential equations (10.3), (10.5) that we may take 
a = —1,b= 1. If, then, we fix k so that P,(z) = 1, we must have 


1 
[ dg = k, or k = 3. 
=], 
Returning now to the familiar notation, and writing n, m for Q, Q’, 
we see that (10.8) is the familiar result 


, 28 nn 
| P, (2) Pin(2) dz = on re) ; 
—1 


We have thus derived the orthogonal property of the Legendre 
polynomials directly from the differential and difference equations. 
+ This should be compared with the work of Ferrar, Proc. Edinburgh Math. 
Soc. (2) 2 (1930), 71. He establishes a connexion between the orthogonality 
and the differential equation, but scarcely such a simple one as the present, 


since he does not keep the differential and difference equations together as 
a conjugate pair. 
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11. Parabolic cylinder functions D,(z). These+ satisfy Weber’s 
differential equation 
(Fat nti) Dye) == @ (11.1) 
and the recurrence relation 
Dj (z)—42D,(z)+D,,4:(z) = 0. (11.2) 
Translated into operational symbols, with g, Q for z, n respectively, 
these equations become 
(p?—39°+3+)x = 9, (11.3) 
(p+4q—eP)x = 0. (11.4) 
The adjoint equations are, then, 
(p?@—34¢°+3+Q)x* = 0, (11.5) 
(p—3q+e-")x* = 0. (11.6) 
It is evident from the self-adjointness of (11.3) that we may take 


x* = y(Q)x, where »(Q) is a function of Q alone. It is then not 
difficult to show that y(Q) must satisfy the recurrence relation 


Qy(@) = W1Q—}). 
Hence we may take y(Q) = 1/T(Q+1). 
The general theorem, Theorem VII, indicates therefore that D,(z) 
and (1/n!)D,(z) are orthogonal functions. So we must have 
b 


| D,(2)D,,(z) dz = kn!8 (11.7) 


m nm 


provided suitable ere a, b can be found. But since we have 
D,(z) = e-**, we may compare (11.7) with 
[ et dz = V2z. 
Therefore the precise form of (11.7) is the well-known result 
[ Dy(e)Dyle) dz = J2rn! 


io 8) 


12. Bessel coefficients /,(z). These satisfy Bessel’s differential 


equation (a? 


ld n 
matsz (2 —5))y (2) = 0 (12.1) 


+ Whittaker and Watson, Modern Analysis, 347. (1927.) 
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and the recurrence relation 


J, -1(2) +S aa(z) = (2n/z)J,,(z). (12.2) 
These we shall write in the equivalent forms, with g, Q for z, n 
respectively, (qpqp+q2—@Q?)x = 0, (12.3) 
{29-41 Q-l(eP + e-P)}y = 0. (12.4) 

The adjoint equations are 
(pgpq+q?— Q?)x* = 0, (12.5) 
{2qg-2— (eP +e-P)Q}y* = 0. (12.6) 


It is immediately clear that (12.5) is satisfied by 
x*(q, Q) = *0(Q)x(q, ®), 
where o(Q) is a function of Q alone, since this transforms (12.5) into 
(12.3). This will further give a solution of (12.6) if 
(2q-40(Q)—(Q+1)40(Q+ Ne” —(Q—1)40(Q—le-P}x = 0 
is equivalent to (12.4). This is the case provided 


o(Q—1)_ o(Q) _ o(Q+1) 


Q—1 @ Q+1 
So we may take o(Q) = Q. 

Theorem VIII then shows that x, x* are orthogonal and normalized 
provided suitable eigenvalues of @ exist for a suitable range of 
integration. Thus it suggests the existence of a relation of the form, 
when we revert to the more usual notation, 

b 


[ (m/z, (z)In(2) dz = KB inn. (12.7) 


Actually there does not seem to be any way of getting the precise 

relation without recourse to the ordinary theory. This gives it ast} 
dz 8 

| Jameel ensalt) = = ee. (12.8) 


ss 2 2(2n+1) 


ie) 


ny 


0 
So the eigenvalues of Q appropriate to the problem are the odd 
integral values. 

This is a case where our theory fails to give the details of the 
orthogonality relation. Nevertheless it does suggest the existence 
and general form of a relation which the detailed theory must seek 
to establish. 


+ Watson, Theory of Bessel Functions, 404. (1922.) 














ON THE SUMMABILITY OF DERIVED SERIES 
OF THE FOURIER-LEBESGUE TYPE 
By A. H. SMITH (Providence, RI.) 
[Received 5 September 1932] 

1. Introduction.; In two recent paperst L. 8S. Bosanquet and 
E. H. Linfoot introduced a method of summability which is a 
generalization of Riesz’s arithmetic means of order « and is defined as 
follows. The series } a, is said to be summable (a«,8) to S, where 
either « > 0, or a = 0, B > 0, if 

y\% 
S Et —- ] log P( 
i n 1 
for C sufficiently large§, where B = (log C)®. When a = 0, B > 0, 
the method has been called zero order summability. 

Throughout this paper the function f(a) will be understood to be 
periodic and of period 27. We shall assume that it is integrable in 
the sense of Lebesgue, but in the latter part of the discussion (§§ 4-6) 
restrict it to be of bounded variation in (0,27). Its rth derived 
Fourier series is defined to be 

bi i'[v"(a, cos va +6, sin va) | (r even) (1.02) 
v=l1 
and > —i*[y"(—a, sin vx +5, cos v2) | (r odd), (1.03) 
v=1 
where a, and 6, are the Fourier coefficients. In particular the first 


derived Fourier series and its conjugate series are given by 


y [v(—a, sin va +-b, cos vx) | (1.04) 
y=] 


v 
2) 


and > [v(a, cos va +6, sin va) | (1.05) 
respectively. = 
If 
f(x | t) + f(a—t) ey f(x)? ff? 2)(ar)t" 2 
ST eT ee 
f(a) wo, (a, tt 
(2r)! © = (2r)! 


+ The author is indebted to Professor J. D. Tamarkin for many helpful 
suggestions. t Bosanquet and Linfoot (1), (2). 
§ Bosanquet and Linfoot (2). They have shown that it is equivalent to 


, (1.06) 


> 


say ‘for every C> 1’. 








94 A. H. SMITH 
where f(x), f(x),..., f(x) are independent of t, and lim w, (zx, t) = 0, 
t+0 


then f(x) is defined to be the generalized derivative of order 2r of 
f(x) at the point x.¢ If 
f(a+t)—f(a— t) - f (aye 4 Faye . f(x) 
2 1! 3 °”™ (2r—1)! 
S haed { D(a)t* +1 w(x, t)t?” +1 
(2r-+1)! (2r-+1)! 
where f(x), f(x),..., f@*?(x) are independent of ¢, and 


lim w.(x,t) = 0, 
t+0 


then f'"+)(x) is defined to be the generalized derivative of order 
2r+1 of f(x) at the point x. 

The Cesaro summability of the Fourier series, its conjugate series, 
and the series defined by (1.02), (1.03), (1.04), and (1.05) has been 
investigated by various writers. The aim of this paper is to investi- 
gate the summability, according to the method defined by Bosanquet 
and Linfoot, of the series (1.02), (1.03), (1.04), and (1.05). The main 


results are as follows. 


(1.07) 


THEOREM 3.1. The rth derived series of the Fourier series of a func- 
tion f(x), integrable-L, is summable (x,B), where « =r, B > 1, to the 
rth generalized derivative whenever it exists. 

THEOREM 5.1. The derived series of the Fourier series of a function 
f(x) of bownded variation is summable (0,8) almost everywhere, for 
B > 1, tof'(z). 

THEOREM 6.2. The series conjugate to the derived series of the 
Fourier series of a function f(x) of bounded variation is summable 


(0,8), for B > 1, to » 
P= tm —* [ s ) at 


n->9 7 


” 
almost everywhere, where x(t) is defined by (4.01). 
2. We shall first define the following functions: 
“Hy 
H,..9(1—u) = Bu!(1—u)*-1 log-8 | — 
tap(l—u) = Bul(1—u)*-og-*( 
forl >0, «a >1, B > 0, where B = (log C)f, 


1 
Qi,0,8(t) = | A, .g(1—u)cos tu du, 


0 
+ de la Vallée Poussin (3), 214. 
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1 
Qro,p(t) = [ H,,.,g(1—u)sin tu du, (2.03) 


0 


1 
= Mapl)+iQ,ap() = [ Hyapll—ude™ du, (2.04) 
0 


A,,g(”, t) = ‘f Hoa(} —*}cosvt dv, (2.05) 
0 


and A‘*}(n,t) will denote the kth derivative of A, 4(n,t) with respect 
to t. We shall understand throughout that K, either with or without 
subscripts, denotes a positive numerical constant. 


Lemma 2.1. For l>0,a>1, B >0, the functions Q,, g(t) and 
Q,..g(t) are bounded in (0,00), and for large values of t 


x 


M baa om 84. 114-0 (2.06) 
wel pa os T \t" loge)’ ies 


where C),. 18 a constant. 

Proof. The first part of the lemma follows immediately as a con- 
sequence of the definitions of Q(t) and Q,, (¢). The second part 
is obtained by using Cauchy’s theorem to express M in terms of 
contour integrals, taken round the rectangle of which the vertices are 
0, 1, 1+7R, and 7R, letting R > oo to derive the expression 

== ft H,,.g(1—u)e™ du — J ar (l—uw)e” du 
0 

and investigating these integrals pide § 

LEMMA 2.2. Aig, t) is bounded in (0,00) for fixed n, where k > 0, 

- 1, B > 0; and for large values of t 

1 

t**+1 logBnt 


1 % = 
= O(a) (« > k+2), 


where k > 0 and B > 0. 
Proof. The boundedness of A‘? ,(, t) follows from its definition. It 


M(n,t) = o( ) (a = k+1), 


+t Cf. Bosanquet and Linfoot (2). 
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is easily verified that 


iknk 


1 
AER(Ns t)= Qie,o,p(Mt) (k even), 


7 


gkt+ink +1 


Q,. ,g(nt) (k odd). 


The lemma now follows from (2.06). 


LemMA 2.3. Whena > 1, B > 0, 
A it = Bix(1—Y)~ 
2 (27)(n, t)cos vt dt = Bi?*| 1— 
| ap )cos vt d ( “} 
0 
- 0 (v : 
and 


2 | Ae) (n, t)sin vt dt 


0 F = 
Bie+(1—*) log | 
n 


= 0 (v >n). 
Proof. The statement follows from the Fourier repeated-integral 
formula, since, for k > 0, a > 1, B > 0, 


»\a-—1 Y 
T.(n,v) = (1 — log A( . ) 
n 1—v/n 


is of bounded variation and continuous for 0 < v < n, where n is 
fixed, and converges to zero as v > n. 


3. We shall now consider the rth derived series of the Fourier 
series of a function f(x), integrable-L, at a point where the rth 
generalized derivative exists. It follows from Lemma 2.2 that for 
fixed n, Apts t) is absolutely integrable and of bounded variation 
in (0,00), when k > 1,«=k+1, 8 >0. Hence, by a theorem due 
to W. H. Young} and Lemma 2.3, we obtain 


2 


| [fle+t) +fla—t) eR (n, t) dt 


y\ a—1 f y 
4 log-* ¢ }»™(a, cos var +b, sin xz) | (3.01) 


mn 1—v/n 


+ Hobson (4), 583. 
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where r > 1, a = 2r+1, B > 0, and 


ao 


— | [fle+t)—fle—rey(n, 1) at 


0 
' pot C ‘ 
—42r+2 neh —fi_ _™~ _\ar+l¢__ 
i p [a(2 “) log (<,) (—a, sin vx +b, cosy), 
(3.02) 


where r > 0, a = 2r+2, 8 >0. Thus we have the following lemma: 


v<n 


LEMMA 3.1. The eapression 


L,,(f,x) = (—1y [ f(x+t)+(—1)fla—t) AY). p(n, 4) dt (3.03) 
0 
is the |n|th} mean of order (r,B > 0) of the rth derived series of the 
Fourier series corresponding to f(x). 


THEOREM 3.1. The rth derived series of the Fourier series of a func- 
tion f(x), integrable-L, is summable (x, 8), where « =r, B > 1, to the 
rth generalized derivative whenever it exists. 

Proof. Investigate first the derived series of even order. Then 


(3.03) becomes 


alf2) = f [e+t)+fle— AE, plr,t) dt (3.04) 


Consider only points where f(z) possesses a generalized derivative of 
order 2r, and thus (1.06) is satisfied. Choose « positive, then h 
positive (but fixed) so that 

lw,(a,t)i<<« (0<t<h), (3.05) 
and then so that nh > e. By Lemma 2.2 


« 


Ly = | [fle+t)-+fee—t) Age, (nt) de 


h 
re 


= O(lou-Pnk | fee) it), 
h 


f2r+1 


thus lim L, = 0 (rf >1,B > 0). 


no 
Next let us investigate 
h 
L, = | [f(a+t)+f(e—t) PAS? p(n, t) dt. 
0 
t+ Where [n] is the largest integer < n. 
H 
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By means of (1.06) set 
L, = 1,4-L,4+L,;, 
where : h ) 
_ > Fora) Qvx(2r) /, 
Ls => on [ AZr) p(n, t) dt 
v=0 0 
fer) r 
2f°”(x) mre 
L= nr) | PAST") 1 p(m, t) dt b 
0 
h 
- 2 / 2r)(2r) 
L, = @n) | wo, (x, t)PrAS") 1 p(n, t) dt 
0 J 


After integrating by parts each of the r—1 terms of ZL, and 
observing that the integrated terms vanish as noo, for B > 0 
(Lemma 2.2), we have 


lim L, = 0 (8 > 0). 


Consider L,. After integrating by parts 2r times and observing 
that the integrated terms vanish at both limits as above, we see that 
the discussion of L, reduces to that of 


(r(x 
=2- 2 fa (n,t) d 


Replace the variable ¢ by t/n and evaluate the integral by Lemma 2.3. 
Then for 8 > 0, lim L, = f®(c). (3.11) 


n> 


It remains for us to investigate L,. 


(3.09) 


(3.10) 


On account of (3.05) 


t?r t\| 
(2r)[, 
ar Xe (5) 


Divide the interval (0,nh) into (0,/) and (l,nh) and investigate the 
two integrals independently. 
¢2r 


—_ (2r) — af 
L, =| ra ) dt<K | iogs = 0(1) 
l 


uniformly in », when a = es and 8 > 1. 


L, val Ae B | 
= ner 1 


dt (a = 2r+1). 








nh 
|L;| < Ke | 
0 


nh 


e 





Finally dt 
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is bounded in (0,00), for all nm and a > 1, B > 0, since 


ATA (Ms t/n) = nao, «Bl ). 
We have now shown that for all values of m and for 8 > 1, 


|L;| < «Ky < «, (3.12) 


where e, can be made arbitrarily small with h. Thus, summing up 
the results of (3.06), (3.07), (3.09), (3.11), and (3.12), we see that 


for B > 1, lim Ly »(f,2) = f(z). (3.13) 


Combining (3.13) with Lemma 3.1 completes the proof for the case 
of the even derived series. The proof for the odd derived series is 
similar. 


4. Throughout the remainder of this paper we shall consider f(x) 
to be of bounded variation in (0,27). Integrals of the Riemann- 
Stieltjes type will be employed. 

We shall define the following functions: 


b(t) = 3[f(x+t)—f(a—t) —2¢f’(a)] | (t) = [f(@+t)+f(e—t)—2f(@)] 
A ‘ (4.01) 


(u) = | (dg(e)| ¥(u) = | ayo) (4.02) 
9p(t) = Qo,1,8(4) p(t) — Qor(?) (4.03) 
A, (x) = v(—a, sin vx +5, cos v2) B,(x) = v(a, cos va +-b, sin va) 

(4.04) 





Lemma 4.1. For B > 0, gg(t) and ga(t) are continuous functions of 
t in the interval (0,00), and, for B > 1, they are of bounded variation 
in the same interval. 

Proof. Continuity in (0,00) follows from the fact that, for B > 0, 
both log-f{C/(1—u)}costu and log-P{C/(1—w)}sin tu are continuous 
functions of (t,) in the region defined by 


O0<u<l, 0<t<o, (4.05) 


and that gg(t) and g(t) converge to zero as ¢ > 0. Bosanquet and 
Linfoot} have demonstrated that for 8 > 1, gg(t) is of bounded varia- 
tion in (0,00). The proof for the case of g,(t) is similar. 


t+ bosanquet and Linfoot (1). 


np 723353 A 
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Lemma 4.2. For B > 0 





[ yp(theos tx dt = Blog-#(5<) (0<2< 1), 


0 


A iv 


= 0 (x > 1); 
and a 
2 ; C 
_ | g(t)sin ta dt = Blog-*( 0<2=— 4) 
7 1—2x 
0 
= 0 (2 = G2 > Je. 


Proof. The lemma follows from the Fourier repeated-integral 
formula since for 8 > 0, Blog-*{C/(1—x)} is of bounded variation 
and continuous in 0 < 2 < 1, and tends to 0as2— 1. 


5. We shall now consider the derived series of the Fourier series 
of a function f(#) of bounded variation in (0, 27) at a point where 


(i) f(x) and f’(x) exist, (5.01) 
(ii) | |\dd(t)| = o(1), as u> 0. (5.02) 
0 


It is known that these two conditions are satisfied almost every- 
where.} We shall first demonstrate the following lemma: 

Lemma 5.1. If f(a) is of bounded variation, and f(x) and f'(x) exist 
at the point x, then for B > 1, 





L 


C ; ee 
> Blog ( . )Ase)-f (x) = lim - | gp(nt)n df(t), (5.03) 
- 1—v/n Loo 7 
van 0 
where A,(x) is the general term of the derived series and the integral is 
taken in the Riemann-Stieltjes sense. 
Proof. First we shall show that for 8 > 1, 


L 





im | gg(nt) dl 3{f(~+t)—f(~—t)}] 
L 
= —lim | 3[fa+t)—flw—t)] dgg(nt). (5.04) 


+ Consult Plessner (6) in the case of (5.02). 











ON THE SUMMABILITY OF SERIES 


Now 


[ gplnt) aa fee+1)—fle—t)}] = [gplne) fe) fee 


0 


L 
— f fle+t)—fle—t)] dgg(nt). (5.05) 


Pass to the limit as L > 00. Because of Lemma 4.1, the integral 
[ ALfe+9—fe—t)] dgp(nt) 
0 
exists} for 8 > 1. The integrated term of (5.05) vanishes at the 
lower limit with [ f(x+t)—f(x—t)], and at the upper limit with 
ge(nt), since it is O((nt)log-Pnt) by Lemma 2.1. Thus the second 
integral of (5.04) exists and the relation (5.04) is satisfied. Similarly, 
using the fact that sin vt is an absolutely continuous function of ¢ in 
(0, L), we obtain 
L L 
lim | gg(nt)A,(x)cosvt dt = —lim | A,(x wo dgg(nt). (5.06) 
[> I> 
0 0 
Next 
L 


L 
: | galnt)n dp{t) = =n | gant) d[3{ fla+t)—fla—t)}]— 


0 


L 
= nf (x x) | gp(nt) dt. (5.07) 


Pass to the limit as L-> oo. By Lemma 4.2 the last term becomes 
—f'(x). Moreover, as L +o, the first term on the right of (5.07) 
exists for 8 > 1, and thus under this assumption 
. i 
aglntym det) = =n | aplnt) A fle+0)—fe—O}]-F'@). (6.08) 
0 0 
Now since f(x) is of bounded variation in (0, 27), 
— A, (a)sin vt 


Vv 
v=1 


converges boundedly to }{ f(~+t)—f(a—t)]. Combining (5.04) and 
+ Hobson (5), 545. 
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(5.06), integrating term by term,{ and then substituting from Lemma 
4.2, we obtain 


=n [ ga(me) d[3{f(x+t)—f(x—t)}| = > =n | gp(nt)A,(x)cos vt dt 


= > Blog-4( :- ),(0) (5.09) 
ven 1—v/n 


From (5.08) and (5.09) the lemma follows. 


THEOREM 5.1. The derived series of the Fourier series of a function 
f(x) of bounded variation is summable (0,8) almost everywhere, for 
B > 1, to f'(x). 

Proof. Assume that at the point 2 the two conditions (5.01) and 
(5.02) are satisfied. Choose ¢ arbitrarily small, then A so that 


u 


- } db(t)|<<« (0<u<A). (5.10) 


Next choosen so that nA > e, and define J,, J,, and J, as follows: 





eln > 
9 
L=- | gg(nt)n ddp(t) 
7 
0 
aan 
I, =- | gp(nt)n dd(t) >. (5.11) 
7 
L 
2 
I, = - | gp(nt)n d¢d(t) 
7 
A 
Since gg(nt) is bounded (Lemma 2.1), then by (5.02) 
lim J, = 0(1). (5.12) 


Next let us investigate J,. 
a % 
al < Ke | opty (M40: 
ejn 
Integrating by parts and observing that the integrated term vanishes 
for B > 1 as n> ©, we see that the discussion of J, reduces to that 


+ Young (7). 











ON THE SUMMABILITY OF SERIES 
of the following expression: 


IK joc (t)tlog-Pnt dt — BK, Jo O(t yHlog-#-In at (5.13) 


In the case of the first integral of (5.13) we may write, by (5.10), 


A A 
-K, | O(t)t-*log-Pnt dt| < Kye tlog-Pnt dt 
ejn 
and for 8 > 1, this term is arbitrarily small with ¢«, since nA > e. 
Also 
|—BK, fom )t-*log-F—-1nt a - Kype [ tlog-F-1nt dt 
ejn e/ n 
and similarly, for 8 > 1, this term is arbitrarily small with e«. 


Finally consider J,. 
L 


[Is] < Ky | tlog-Pnt |dgi(0) 
A 
On account of the periodicity of f(x), if g and s be chosen so that 
2(q—l)xw# < A < 2qn < 2As—1)z < L < 2Qsz, 


Qn s—1 


1 1 1 
then |Is| <K f [ aa \ddi(t) ai log?nA * On >) “ioe 


wo 


The series pA a “oe converges, for 8 > 1, to a sum which tends 
to zero as n> 00. Thus for B > 1, 
re 
lim = [ gp(nem dd(t) > 0, as n > 0. (5.14) 
Lox 7 
In view of the fact that the conditions (5.01) and (5.02) are satisfied 


almost everywhere, the theorem follows from (5.14) and Lemma 5.1. 


6. We shall next investigate the summability of the series con- 
jugate to the derived series of the Fourier series of a function f(z), 
of bounded variation in (0,27), at a point where the following con- 
ditions are satisfied: 


(i) a |deb(t)| >0, as u>0, (6.01) 


0 
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) lim— fe. t exists. (6.02) 


70 


We state first the following ieee: 


L 
lim [ gg(nt) d[3{fle+t)+fe—t}] 
Lox’ 
0 
L 
= —lim | al s[ f(w+-t)+f(~—t)] dgg(nt), (6.03) 
Lox” 
L ‘ L 
, E : ; B,(«) a 3 
lim | Jp(nt)B,(x)sin vt dt = lim | —*—cosvt dgp(nt). (6.04) 
I-x Lo Vv 


0 0 
The existence of the Riemann-Stieltjes integrals and the proofs of 
the relations (6.03) and (6.04) are similar to those given for (5.04) 
and (5.06). 





Lemma 6.1. If f(x) is a function of bounded variation, then for B > 1, 

L 
> Blog- ‘(< iq) B (x) = lim —2 | gglnt di(t), (6.05) 
v<n — o 3 

where B,(x) is the general term of the series conjugate to the derived 

series and the integral is taken in the Riemann-Stieltjes sense. 


Proof. We have 


is 9) is) 


2 | Gplmt) dM fla+t)+fw—t}] = | Gp(nt) adie) 


0 0 
The remainder of the proof of the lemma is similar to that of 
Lemma 5.1. 


THEOREM 6.1. The series conjugate to the derived series is summable 


(0,8), for B > 1, to 
= lim — 1 dt, 


n>0 


whenever this limit exists, provided es (6.01) is satisfied and b(t) = o(t) 
ast—>0. 
Proof. It is sufficient to show that under these conditions 
L 
J =lim—” | Gp(nt) df(t) > T, as n> 00. (6.06) 


I-7n 7 
0 
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Choose ¢ arbitrarily small, then A so that 


1 J \d(t)|<e (0<u<A). (6.07) 


Next choose n so that nA > e, and define J,, Jy J, as follows: 


eln 


i - =” [ ggint dib(t) 
7 
0 





In the interval (e/n, L), 
Gp(nt) = (1/nt)+hg(nt), 
where by (4.03) and Lemma 2.1 
hg(nt) = Of(nt)log-Fnit}. 
Set J, = J+J and J, = J+J." 
where r 
t= [cm (nt) dip(t) 
7 
oa (6.09) 


Ji = —* | hg(nt) dip(t) 





with similar definitions for J* and J**. As in Theorem 5.1, the 
integrals J,, J**, and J** are o(1) as n> oo, for 8B > 1. Thus for 
n sufficiently large and f > 1, 

L L 
lim —" | gy (nt) dyp(t) = lim —4 fe di(t) + 0(1), (6.10) 


L— Lo 7 
0 ejn 


where if dip(t) = ~(2)4(3)—= | t-2yp(t) dt 


Let n > oo, then, since (t) = o(t) as t> 0, 
tim —> [ 1 dip(t) = - a if t—1 dy(t) = lim at h(t) 


n> oO n->0 
e/n v v) 


provided the limits exist. Thus the theorem is proved. 
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Let us now consider the restrictions of Theorem 6.1. It is known 
that (6.01) is satisfied almost everywhere.f Moreover, since ‘(é) is 
the total variation of #(¢) in the interval (0,¢), it follows that 
%(t) = o(t) as t> 0 almost everywhere. Plessner{ has established 


the existence of “ 
lim [ WO) ay 


no J sin?(3t) 
7 
except at a set of points x of measure zero. By writing 


1f we) =. it) 
i] sin?(4t) - ak | esa $t) Jas [te 


: 2 
the existence of lim — — 
n->0 7 

i] 

almost everywhere is demonstrated. Hence we have the following 


F spit) 
‘2 


theorem. 


THEOREM 6.2. The series conjugate to the derived series of the 
Fourier series of a function of bounded variation is summable (0, 8) 
for B > 1, to 

T =lim— f ae 


n->0 7 
n 


almost everywhere, where x(t) is defined by (4.01). 
+ Plessner (6). 
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THE FORMAL REDUCTION OF SYSTEMS OF 
LINEAR DIFFERENTIAL EQUATIONS 


By E. G. C. POOLE (Ozford) 
[Received 16 October 1932] 


Introduction 


In his treatment of the classical problem of systems of ordinary 
linear differential equations with constant coefficients, Jordan gives 
a remarkable elementary process* for exhibiting the invariant factors 
of the system and reducing it to an equivalent diagonal system of 
special form. So far as I know, it does not appear to have been 
noticed that the same method is applicable with slight modifications 
to systems with variable coefficients. 
Consider the system 


U; = > Fis y; == § (a = 1, 2,...,”), (1) 
f=1 
where (F;;) are operators of the type 
q™ qm-1 
F= o(%) 5 on + 0,(#) Tana + ...+4,,(2), (2) 


of maximum order m, say. We shall show that (1) is equivalent to 
a diagonal system 

H, w, = 0, Hiw,=6, .., Hw, = 0, (3) 
where some of the operators (H,) may reduce to unity. Further, if 
the equations (3) are arranged in the correct order, they have the 
property (analogous to that of invariant factors) that every solution 
of Hw = 0 is a solution of H;,,,w = 0, and every integrating factor 
of H;.w = 0 is an integrating factor of H,.,w = 0. 


Inner and outer factors of operators 
If two equations (F,y = 0, F,y = 0) have common solutions, we 
may perform Brassine’s process} (analogous to the extraction of the 
highest common factor of two polynomials) and find (K,y = F; Gy, 
F,y = F; Gy), where Gy = 0 is the equation satisfied by the common 
* Jordan, Cours d’analyse, iii. 174-9. 


+ See, for instance, Schlesinger, Theorie der linearen Differentialgleichungen, 
i. 42-5, or Picard, Traité d’analyse, iii. 554 (ed. 1908). 
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solutions. We may say that the operator G is the highest inner 
common factor of (F,, F,). Suppose now that we require the integrat- 
ing factors common to the two equations; the classica] solution is to 
construct the two adjoint equations and apply Brassine’s process 
to them. We may, however, proceed more directly, if we observe 
that any common integrating factor of (F,y = 0, F,y = 0) is also 
one of F, P,y+F,P,y = 0, where (P,, P,) are any operators. By 
an obvious counterpart of Brassine’s process, in which the order of 
products of operators is reversed, we reduce the equations to the 
form (HF{ = 0, HF = 0), where Hy = 0 admits every integrating 
factor common to the given equations. We shall then call H the 
highest outer common factor of the operators (F,, F,). 


Reduction by rows and by columns 

The passage between the equivalent systems (1), (3) is effected by 
a sequence of elementary steps of two kinds. In the first place, a pair 
of equations (u; = 0, u; = 0) may be replaced by (u;+Au; = 0, 
u; = 0), where A is any linear operator; for either pair is a conse- 
quence of the other. If we perform a sequence of these steps and 
keep a record of what we have done, we obtain two associated 
operator matrices connecting (wu; = 0) and (v; = 0), such that 


n n 
uu, = 2 Ai v,= 2, Ais %- (4) 
j= j= 


By substituting one set of formulae in the other, we see that the 
operators satisfy identities 
- . 
2 Aw Ars = 84; > Ain Ans = bi, (5) 
where 5;; = 0 (¢ AJ) or 8,; = 1(¢ = 7). This process we call reduction 
by rows of the given system. 

But we can also make changes of variables, using instead of the 
pair (y;, y;) the pair (y;+-Ay;, y;), where A is any linear operator; we 
see that either pair can be explicitly written in terms of the other. 
The sets of variables (y;), (w;) are connected by a sequence of such 
steps, and, on eliminating intermediate terms, we find again two 
associated matrices of operators (B;;), (Bj;) such that we have, sym- 
bolically, (y = Bw, w = B’y), as a shorthand for relations-of the 
type (4) and (5). 
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Generalized Jordan process 

First Stage. 

For simplicity of exposition, we divide the reduction into two 
stages; we obtain first an equivalent system 

G,z, = 0, G.4,=6, .., Gz, = 6, (6) 
which does not in general satisfy the supplementary conditions im- 
posed on (3); it is then sufficient to consider the equations (6) pair by 
pair and show how they may be further reduced to the required form. 

We shall suppose that a system has been ‘simplified’ if it has been 
shown to depend on one with fewer equations; it does not follow that 
the ‘simpler’ system is necessarily easier to solve. Further, of two 
systems having the same number of equations and variables, we 
define that one as the ‘simpler’ which has an operator (not identically 
zero) of lower order than any operator in the matrix of the other. 

(i) If any operator of the matrix (F;;) is of order zero, but does 
not vanish identically, we can express the accompanying variable 
explicitly in terms of the rest and eliminate it from the system; the 
new system has fewer equations than the old and is, by definition, 
simpler. For example, to fix the ideas, let F,, be unity; then y, is 
known in terms of (Ys, ¥3,...,¥,). If we write 

Wy = IW t+Foyot--t+Fin yn (7) 
one of our new variables is w, and one of our new equations is 
H,w, =w,= 0. In practice, we merely eliminate y,, keeping a 
record of its expression in terms of the variables of the simplified 
system. 

(ii) If no operator (F;;) is of order zero, we select the operator of 
lowest order of the matrix, or one of the set of the same minimum 
order. To fix the ideas, let us again suppose that it is #,,. We now 
try to reduce the system by rows to a simpler equivalent one; to do 
this we examine whether F,, is an inner (or right-hand) factor of the 
other operators in the same column F,,, F3;,..., F,;. If it is not, 
suppose that F,, is not divisible: on the right by F,; we have 
F,, = F, Fy,+Ry, say, where R,, is an operator of lower order 
than F,,; we can now replace u, = 0 by —F}, u,+u, = 0, leaving 
the other equations unchanged. The modified system is evidently 
equivalent to the original one, and, by definition, is simpler. We note 
that R,, is not necessarily the operator of lowest order in the new 
system, since the new equation may happen to contain a lower one. 
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If at any stage an operator of order zero appears, we reduce the 
system by eliminating a variable, as in (i) above. 

(iii) If we do not see any possibility of reduction by rows, we 
attempt a reduction by columns; to do this, we examine whether our 
operator of lowest order, F,,, is an outer (or left-hand) factor of each 
operator in the same row, u, = 0. If it is not, suppose that F,, is 
not divisible on the left by F,,; we have F,, = Fy Fj,+Sj., say, 
where S,, is of lower order than F,,; introducing a new set of 
variables (y,+ Ff, Ys, Yo, Yg:---» Yn) equivalent to the old, the system 
is again simplified. 

(iv) We continue to apply processes (i), (ii), (iii) until no further 
reduction appears possible. Suppose that the new system* still has 
nm equations and n variables, and that the operator of lowest order 
is again F,, (which does not imply any loss of generality, as we have 
only to rearrange the order of our equations and variables). By 
hypothesis, /,, must now be an outer factor of every operator of 
the new equation, which we call u, = 0; this can be written accord- 
ingly in the form 

Wy = Fy% = Flt Freyot Fisyst ---+Finyn| = 9, (8) 
and defines a new variable z,, such that the set of variables 
(21, Yos-++> Yn) i8 equivalent to (y,, Yo,...,Y,). We can now eliminate y,. 

But F,, must also be an inner factor of every operator in the first 
column; this column of the matrix is unaltered by the change of 
variables, and its operators are of the form F3, F,,, Fj, Fy,..., Fis Fy 
We can now replace the equations (wu; = 0) by the equivalent set 

u, = 0, — Fi, u,+u, = 0, 
—Fy%+uUg=0, ..., —Fh,u,+u, = 0. (9) 

The variable z, now appears only in the first equation of the new 
system; we have now 

Gz, = Fz, = 0, (10) 
together with a system of n—1 equations in Yp, ¥3,..., Yn- 

By applying the same methods to this residual system, we can 
gradually evolve the diagonal system (6). 

Second Stage. 

We may set apart any of the equations (6) which are of the form 
z; = 0, no further reduction being necessary. Let us consider any 


* For simplicity, the notation (1) is retained schematically for the system 
at this stage of the reduction. 
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other pair of equations, say, 

G,z, = @, G,2z, = 0, (11) 
where G, is of lower order than G, (or of the same order). Unless 
the operators already possess the characteristic property of the 
operators H;, namely that G, is both an inner and an outer factor 
of G,, further simplification is possible. 

In the exceptional case, we shall have two identities of the form 
Gz = G,G,z, Gz = G4, Gz. (12) 
We observe that it does not necessarily follow that G, can be written 
in the form G,G,, as may be readily verified by considering opera- 
tors with constant coefficients. 
For suppose that G, is not an outer factor of G,. We replace (11) by 
the system 
G2,+G,z, = 0, Gz, = 0. (13) 
If G, = G,G+8,, we take new variables (z,+ G}z,., 2.) and write 
G, (2+ Go2.)+ 8,2, = 0, Gz, = 0. (14) 
Since S, is not zero and is of lower order than G,, the system is, by 
definition, simpler than before; and as it is not of diagonal form, we 
can resume the sequence of reductions by the standard methods. By 
considering the process for the extraction of the highest common 
factor, it can be shown that the sum of the orders of the equations 
in two equivalent diagonal systems is the same; this point may be 
left to the reader. 
Again, suppose G, is not an inner factor of G,. We replace the 
variables (21, 2.) by (2, 2.—2,), and the equations become 
G,z, = 0, G.2,+G,2, = 0, (15) 


where 2, = %,—2,. The system (15) is reducible by rows; for if 
G, = G,G,+ R,, we operate with G} on G,z, = 0 and subtract from 
the second equation, obtaining the equivalent system 

G,z, = 0, R,2,+ G2; -= 0. (16) 


Since R, is not zero and is of lower order than G,, so that (16) is, by 
definition, simpler than (15). The reduction of (16) to a diagonal 
form may be completed as before. 

Unless (G,, G.) fulfil the conditions (12), simplification is always 
possible, and in the new pair of equations, one will be of lower and 
one of higher order than either of the old pair. For example, if 
(G,z = 0, G,z = 0) have no common solution, and are of order 
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(n,, %2), we can obtain an identity 
K,G,2+K,G,z =z. (17) 
Also, from the equation of order n,-+-n, satisfied by every solution 
of either G,z = 0 or G,z = 0, we have the identity 
Iz = P,G,z = P, Gz. (18) 
The general solution of Iz =0 is z=2,+2,, where G,z, = 0, 
Gz, = 0. But we can also express (z,, 2.) in terms of z; for we have 
z, = (K,G4,+K,G,)z, = K, G2, 
or z, = K,G,(z,+2.) = K,G,2, (19) 
so that the solution of the system (11) is 
2, = K,G,z, 2, = K,G,z, (20) 
where z is any solution of the single equation Lz = 0. In the com- 
plete process, a second variable w would appear, such that (z,, 2.) 
and (z,w) are equivalent, and we should have the equation w = 0, 
forming with Lz = 0 the reduced system. 

Thus any pair of the equations (6) can be replaced by an equivalent 
pair, having the special property (12) that one is both an inner and 
an outer factor of the other. This is the characteristic property of 
the operators (H;) of (3). So long as any two operators of the system 


(6) do not possess this property, a further reduction is possible. In 
a diagonal system of the type (6) taken at random, it is improbable 
that the operators will have common inner or outer factors, and we 
shall usually find that all the variables z can be made to depend on 


a single variable w. 








THE MOTION OF THE PROBABILISTIC FLUID 
OF PRESENCE OF THE ELECTRON IN 
RELATIVISTIC WAVE MECHANICS 


By U. CRUDELI (Cagliari) 
[Received 26 January 1933] 


THE formula obtained by me* for the probability of the presence of 
the electron, with regard to the relativistic wave mechanics prior to 
Dirac’s modifications, allows us to go farther in the study of the 
behaviour of the same probability. 

The following note derives an equation for the motion of an 
electron in an electromagnetic field from a relativistic wave-equation 
(with regard to the new meaning of the motion of an electron). The 
new equation contains certain terms not appearing in the corre- 
sponding classical equation. 

Having written (under a general form) the de Broglie-Schrédinger 


funceon gb = aerridlh (h = Planck’s constant), (1) 


where a and ¢ are both real functions of the space-time coordinates, 
and having substituted (1) in the fundamental relativistic equation, 
which governs the above function, we obtain the two following 
equations: 

a{_|¢—2(U—Grad ¢) x Grada = 0 


h? , (2) 
—— []a+(U—Grad ¢)?+ mc? = 0 
4n7a 


4 
g2 g2 
wheret (j= Zz — —A,+-=>; 
S Ys ys 


‘S4 


with ¥, = 2, ¥=Y, Y3 = Z, Y, = tt, 
and where we have introduced the two four-dimensional vectors 


ene be Wee Te ty 
OT?" Fe ee 


a (F op op “| 
Grad ¢ = (—, —, —,— 
PY, YY2 Ys Ya 
* Presented at the International Congress of the Mathematicians in Zurich 
(September 1932). See Verhandlungen, ii. 309. 
+ Cf. Gordon, Zeitts. fiir Phys. 40 (1927), 119, equation (9). 
t We use (following the Italian convention) the symbol x as a sign of 
scalar multiplication, and the symbol / as a sign of vector multiplication. 
3695.4 I 
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the constants —e« and m denote respectively the charge and mass, 
and finally V and A denote the scalar potential and vector potential 
relative to an assigned Maxwellian field, in which our electron is 
supposed to be susceptible of free presence. 

The first of (2), by means of Gordon’s four-dimensional vector* 


and of the classical circumstance 


—} > 6 
c ot <1 OX 


leads to the continuity equation 


a 
oP 1 div(pw) 0, 
where p represents for us the density of the probabilistic fluid of 
presence of the electron, while w denotes the (vectorial) velocity 
of the fluid itself. 
Finally, introducing the quantity 
h? a . 
Q or 9.9 —_— c*A,a 
8ar“c-ma \ ot" 

and recalling the second of (2) together with Gordon’s four-dimen- 
sional vector, we obtain our general (exact) formula, mentioned by 


us at the beginning of this paper, 


9 1 92 1 
r= (tenga) 
mc~ c* 


Let us not forget to notice that substituting ka in the place of a 
(with k = constant different from zero) the equations (2) remain 
unaltered; and so we may substitute ka for a in this formula. 

It may be useful to remember that, after having written the com- 
ponents of Gordon’s J four-dimensional vector in the form 


J, = mpw,, J, = mpw,, J, = mpw,, J, = temp, 





we obtaint mpw = a(' A—grad $). (3) 
- 


Let us write, for brevity, 


20\3 y2\_4 
p= m(1+ “(1-3 " 
mc? c 


* Loc. cit., p. 121, equation (14). 

+ See also L. de Broglie, ‘The wave mechanics and the atomic structure 
of matter and of radiation’, in Selected Papers on Wave Mechanics (1929), by 
L. de Broglie and L. Brillouin, p. 134. 
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and keep in mind (3); we shall have 


pw = -A—grad 4, (4) 


the meaning of the symbol grad (written with small initial) being 
evident; and then 





d __€0A od eae ie 
ou) = Gerad | 7) + (wx rad) (ew), (5) 


indicating by (w x grad)(uw) the vector, which has 
Ww X grad(uw,), Ww X grad(yw,), Ww X grad(yuw,) 
for its components. 
Let us notice also that the second of (2) leads to 
(«r+ z) = = m*ct+ (eA—c grad ¢)?+ 2mc?Q. (6) 
D 
From which, taking the gradient, we have 
2 (« i+ “4 grad («1 r+ 4 = c* grad(uw x ww)-+ 2mce* grad Q. 
Cc C 


But 
grad(uw xX uw) = 2(uw x grad)(uw)-+ 2uw /\ curl(uw) 


- 2u(w Xx grad)(uw)+ 2° uw A curlA; 
so that , 


(< V+ “f) grad (eV + 4) 
ot ot 


= c*u(w >» Pee A curlA+me? grad Q. (7) 


From (6) re with (4), we have also 


retoon fees} 


where the radical is assumed by us to have its arithmetical value. 

Now, introducing the hypothesis of the smallness of the ratios w/c? 
and 2Q/mc*, we shall be able to substitute, in (8), for the second 
member the expression 


on (28 


and then replace this by pc’. 
12 
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Therefore, instead of (7), we shall write 


grad (2) = —egrad V+(w x grad)(uw)+ : w (A curlA+ ™ rad Q. 
be 


ot 
Finally, substituting in (5), we obtain 
¢ uw) = e(grad “ee ne —*w/curlA— ™ erad Q 
dt ca} c mm 
or 4 (uw) = ~<(B-+ow /\ H) — "grad Q, (9) 
dt c bp 


with evident meaning of the symbols. This is the desired equation. 
As we can see (let us say in passing) the corresponding classical 


2\_1 
¢ mw (1) 2| _ —<(E+2w AH) 
( } c 


C2 


equation 


which differs from (9), might be obtained merely by taking no 
account of (m/u)grad Q in the second member of (9) and by replacing 
in the first member the quantity » by 


2\ 1 
w*\—s 

m{i1——] °. 
C2 














ON SOME CONGRUENCES ASSOCIATED WITH A 
SUB-VARIETY IN A J, 


By H. A. HAYDEN (Cardiff) 


[Received 9 July 1932; in revised form 14 November 1932] 


1. Introductiont 

Iv is well known that the notions of curvatures and normals of a 
curve can be immediately generalized to relate to a congruence of 
vectors defined along a curve. Thus, if (* is a congruence of unit 


vectors along a curve C in a Riemannian space V,, we can define 

n—1 unit congruences, (%, [%,..., f*, and n—1 scalar functions k,..., k , 
2 3 n 1 n—-1 

as follows:{ kf* is the component of D{*/Ds perpendicular to the 

rr+l r 
r-space determined by the vectors ¢*%,..., (*, where (* = (* and D/Ds 
1 r 1 
denotes absolute differentiation along C, and the k’s are all positive. 
The congruences (%,..., (* are mutually orthogonal and satisfy the 
1 n 
equations 


Doi 


r — k(x _k fo (r = 1,....n—1; k = k = 0). (1) 


Ds rrti r-1r-l 0 n 


(When ¢* is tangential to the curve, these are the Frenet formulae 
and k = x.) 
r r 


We shall call the congruence (* the ‘rth normal congruence of the 
r+l 
congruence (”’ and k the ‘rth curvature§ of the congruence {*’. 


McConnell (6) shows how the curvatures may be expressed in 
terms of ‘rates of turning’ of the normals. 

If C lies in a sub-space V,, of V,, and if the congruence (* is 
there are also defined, in a similar manner, m—1 


tangential to V,, 


+ The notation in §§ 1, 2 is that of (3) and (4). A, denotes an affine n- 
dimensional space, V, a Riemannian space, E,, a Euclidean space. 
+ This simple and direct method of generating successive normals is given 
by McConnell (6). 
§ Eisenhart (2) uses the terms ‘associate directions of order r’ and ‘associate 
curvatures of order r’ (extending Bianchi’s terminology) ; Struik (9) calls k the 
, 


‘rth curvature of the curve C relative to the congruence (*’. There appears 
to be no ambiguity in the terminology used in the text. 
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‘normal congruences’, ¢’¢,..., ¢’4, and m—1 ‘curvatures’, k’,..., k’, 
> 


m 1 


satisfying the relations 


relative to V,,, 
Dei 
r =f (*—Z [% = m—1;k'’ =k’ = 0), 
0 


Ds rrtl r-1r—1 m 
J .  — baat ? ; 
where (’* = (’* = @ C*. 
1 


We shall define a congruence along a curve in V,, as being ‘con- 
jugate of order p’, if its first p normal congruences relative to V, are 
all tangential to V,,; p is essentially not greater than m—1. This 
extension of the term ‘conjugate’ is justified in § 2. 

An asymptotic line of order p is, therefore, the envelope of a con- 
jugate congruence of order p, and the object of the earlier part of 
this paper is to obtain some properties of conjugate congruences 
corresponding to those of asymptotic lines given in (3). In particular, 
it is shown that, for a V,_, in V,, the Beltrami-Enneper formula is 
only a special case of a simple formula applicable to any curve 
in V,_,. 

In the later part (§§ 3, 4) we consider the more general case in 
which V,, (or the portion of it contiguous to the curve) is replaced by 


a non-holonomic variety V7", an important feature of the generaliza- 
tion being that V’" is not restricted to be tangential to the curve. 


2. Congruences tangential to a J, 

By arguments similar to those in (3), 25-7, it is easily shownt 
that each of the sets of conditions [A;p], [B;p], [C;p] below is 
necessary and sufficient in order that the congruence {* be conjugate 
of order p: 


[A P| co — C/19.% 
1 


r+1 r 


[ B; p] k= 


i.e. the first p normal congruences and curvatures relative to V, are 
identical with those relative to V,,; 
[C; p] 0;;°E"C'I = 0 (r = l,...,p; a = l....,m), 
r 
where €” is the tangent vector to the curve, regarded as a vector in 
V,,, and Q;;* is the Eulerian curvature tensor of J, in V,,. 
For a V,_, in V, the conditions [C;1] for a conjugate congruence 


+ The reader will find no difficulty in supplying the proofs. 
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of the first order reduce to the single condition 
Wi; . Yi = 0, 
where w;; is the second fundamental tensor of V,,_,; so that the 
directions of ¢’' and £”' are conjugate in the ordinary sense. The notion 
of conjugate congruences of higher orders is therefore an extension 
of the ordinary notion of conjugate directions. 
Before proceeding further, we consider whether there are any 
restrictions on the existence of conjugate congruences.t 
The V,-vectors 0;;°¢", that is, D0;*/Ds, for j = 1,..., m, at a point 
of C, all lie in the first normal spacet of V,, relative to V,, and 
possibly in a space of lower dimensionality. Denoting the dimen- 
sionality of that space by MU, let x* (o = 1,...,M) be any M mutually 
o 


orthogonal unit vectors in that space; then equations [C;p] are 


equivalent to 


w, Ct = 0 (r = 1....,p; 0 = 1...., M), (3) 
r 
where 
o o . D6;% Dyxq 
} ’ wid , U a -a@ o 
oa; = wy, € = OF Xue" = ~ Xa Ds 7 6; Ds , (4) 
Equations (3) are equivalent to 
e refi 
one = Q (r = 0,1,...,p—1; o = 1...., M), 
8 
which are easily seen to be equivalent to 
. r ae 
pid = x @ (r = 0,1,...,p—1; o = 1...., M), 
Ds* 
that is, to o 
C’tw; = 0 (y = 1.,...,9; @ = 1,..., Mt), (5) 


r 
where w, is the (r—1)th normal congruence of w; relative to V,,. 
The number, N (say), of independent equations (5) is equal to the 
number of independent vectors w; (r = 1.,...,p; o = 1.,..., M) at the 
r 


point considered; thus N < m and also N < pM < pm,, where m, 
is the dimensionality of the first normal space of V,,. For the 
existence of a conjugate congruence of order p along a given curve 


+ I am indebted to Mr. A. G. Walker for suggesting the consideration of 
existence-conditions. t Vide (4), 45. 
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C, it is necessary and sufficient that N < m at every point of C. 
Hence, if m, < m/p at every point of C, there certainly exist con- 
jugate congruences of order p along C. 

In particular, if every point of C is an axial point of V,, (in which 
case m, = 1), there exist conjugate congruences of all orders along C. 
Further, since m, < n—m, we have the following theorems: 

I. Conjugate congruences of all orders less than m/(n—m) exist along 
any curve ina V,, in V,. 

Il. If m > 3n, conjugate congruences of the first order exist along 
any curve in a V,, in V,. 

III. Conjugate congruences of all orders exist along any curve in 
a V,_, tn JV... 

Also, since M = 1 in this last case, the conjugate congruence of 
highest order (viz. n—2) along a given curve is unique, in general. 

The mth curvature relative to V, of a conjugate congruence of 
highest order (viz. m—1), if such exist, is given by the relation} 

kh? = OO EE = OG 2 QF, EE'*—(m—1)H 
m 
= mN, =*—(m—1)H, 
=* is the 
normal curvature vector in the direction of the curve C, and H is 
the mean relative curvature of V,, in V, in the direction of the 


where N® is the mean normal curvature vector of V. 


m> 


curve. Hence _—o ae 
k? = mJvcosd—(m—1)H |, (6) 


m 


where J is the mean normal curvature of V,,, v is its normal curvature 


in the direction of the curve, and @ is the angle between the vectors 
=“ and N*, 


For m = n—1, the directions of the vectors E* and N*“ coincide 


in that of the normal to V,_,, and hence, for the conjugate con- 
gruence of order n—2 along any curve in V,_,, 

o2 = (n—1)Jv—(n—2)H; (7) 

n—1 
while, for a ¥, in #3, the second curvature of the conjugate con- 
gruence (of the first order) along any curve is equal to 

¢ rr 
J(2Jv—K’), (8) 

where K’ is the Gaussian curvature of the surface \j. 


+ Again cf. the proof in (3), § 7. 
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The first term on the right-hand side of (6) vanishes if 
057, Qi EE* = 0. (9) 

If J is zero at a point of V,,, every direction in V,, at that point 
satisfies this equation; thus, for example, the first term in (6) vanishes 
for every curve in a minimal sub-space of V,. 

If J is not zero, the directions é" satisfying (9) generate an (m—1)- 
dimensional quadratic cone at each point of V,,. This cone contains 
the asymptotic directions (when such exist+), corresponding to the 
case v = 0; the remaining directions, which exist only if m < n—1 
and m, > 1, correspond to the case 6 = + 4. 

When the curve C is an asymptotic line, the congruence ¢" is 
tangential to C, k is equal to x, and the formulae (6)—(8) reduce to 


m m 


the classical Beltrami-Enneper formula and its generalizations, (3), 
30-1. 


3. Congruences in a VV” in JV, - 
If BX (i = 1,...,m) are m (< n) linearly independent vector fields 
i 
defined at points of a sub-space J; of V,, the aggregate of the vector- 
spaces determined by them constitutes a non-holonomic variety 
Vi in V,. 
A connexion I%!, ‘induced’ in V7" by V,, is given byt 


je 
ti ; (eBj* e B:AQg- 
mi = Bi, +T%, Bj), (10) 


¢ 


where (i) y¥ (c = 1,...,l) are coordinates in Vj, and 0;4 = dat/dy’; 


(ii) Bj* = é B, where é (j = 1,...,m) are the fundamental covariant 


J 
vectors (‘Massvektoren’) in Vj" corresponding to the system§|| B*, 


and Bi, B;* = 8}; (iii) Py, are the coefficients of the connexion of #8 
If V, denotes covariant differentiation ‘secondo R. Lagrange’ (the 
type of differentiation, applied to tensors having several distinct sets 


+ The number of asymptotic directions is «0”"-™, 2", or 0 according as 
‘ > m—1 (Struik (9), p. 110). 

t Schouten (7); Schouten und van Kampen (8) (for A” in A,); Bortolotti 

(1) (for A®, pseudonormal to A,, in A,). 


m 


my, 


§ That is, e, are the vectors whose components in the given system are 
(1, 0,..., 0); (O, 1,..-) 0); ..-3 (0,..., 0, 1) [cf. (1), (8)]. 
So far as Latin indices i, 7, k,... are concerned, contravariance and co- 
variance is with respect to linear transformations of the vectors Be. 
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of indices, which Bortolotti employs to such advantage?), 


v, Bj = 2Bi* 5 ry B:A0; 414 Bi. (11) 


wy 
This may be called the Eulerian curvature tensor of Vj" in V,, and 
will be denoted by Q;;%. 

Since it is only the portion of V,, contiguous to the curve C which 
enters into the considerations of § 2, we shall restrict ourselves, with- 
out loss of generality, to the case of a V7" formed by m linearly 
independent vector fields Be along a curve C in V,. In this case 


there are m? coefficients of the induced connexion, 


ai - ee” . se ae 
i= Bi, s +I i. Bie) 


where &¥ is the unit tangent vector to the curve; and the Eulerian 
curvature tensor of V?" is 

Bx “0 ee 
Q;« = DBT _ IBS" pe BAgu_y Ba. (13) 

Ds ds ' 
We shall write this as Q;*%. 

If (* is a congruence of unit vectors along C in V” and ¢" are its 
components in V%" (viz. B‘,¢%), we obtain two sets of normal con- 
gruences and two sets of curvatures—one set relative to V,, the other 
relative to V’”"—satisfying equations analogous to (1) and (2). 

If we define a congruence in V’" as ‘conjugate of order p’ (p < m) 
if its first p normal congruences relative to V, all lie in V”, then, 
proceeding in a manner analogous to that in § 2 and remembering that 

BY Q;* = 0, (14) 
as is easily verified from (12) and (13), we find the following necessary 
and sufficient conditions for a conjugate congruence of order p: 
[A;p] C% on £85 (y = .,...,9), 

r+l r+1 
[ B; p| k= k’ (r = 1....,9), 
r r 
[C;p] 0;*f"* = @ (yr = I,..:,9). 
_ 
Also, in the case of a conjugate congruence of order m—1, if such 
Oxiey, k? = Q;2Q',, (15) 


m 
where Qi, = DB',/Ds. 


+ It is fundamentally equivalent to the ‘D-Symbolik’ of van der Waerden 
and Schouten [vide e.g. (8)]. 
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There are restrictions on the values of p analogous to those in § 2; 
while Theorems I-III of § 2 remain true if we replace the phrase 
‘along any curve in V,,” by ‘in any V?"’. 
If B3 is the fundamental tensor (‘Einheitsaffinor’) of V7", viz. 
Bz = > Bx By = B;* B',, and if we write Qf = DB3/Ds, it is easily 
c=1 ¢ i 
verified that [C;p] and (15) take the alternative forms:t+ 


[C; p]’ 03 0 = 0 (r = 1...., p), 


r 


k? = 40308. (15)’ 


If m = n—1, Bs = 83—X°Xp 
where y* is the unit vector normal to V?-1, and we find, for the 
conjugate congruence of order n—2, 
_ Dx* Dxa. 


2 — : 
n—1 Ds Ds 


i.e. k is equal to the first curvature of the congruence of normals 


If, when V” is tangentiait to Vj, we define a curve in J); as an 
‘asymptotic line§ of order p of V7’ if its first p normals relative to 
Vv lie in V7", the conditions for an asymptotic line of order p are, 


therefore, "i R- , 
g* _ é ‘az, 


cox, 
“ 


+1 r+1 


: 
Oe £" f= Ove £"e £8 = 0, 
r rT 
where é”¢ are the components of the tangent-vector in the coordinate 


+ To obtain these, we note that 
Q;* = BgQ;F+B FOZ = BPO, 

whence k? = BaQ* a8. 
Also, since B; B* = Bi, 
it follows that Bs; 08 + BP OS = OX, 
and hence that 2B3 OA OF = OF 0%. 

{ i.e. the tangent l-spaces to V, are contained in Vj' if 1 < m, or contain 
Vr itl > m. 

§ Asymptotic lines in a non-holonomic variety (V;*~' in V,) have been 
considered by Hlavaty (5). 
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system y° (in V;), x’ is the rth curvature of the curve relative to V" 


s 
(not relative to \;), and Qg = V, BB. 

For an asymptotic line of order m—1, 

x? = O7°O3!  E%E"4 = 4052054 E"E"4, 
m 
4. Congruences normal to a V" in V, (m < n—1) 

There is an (n—m)-direction at each point of C perpendicular to 
V*, and their aggregate constitutes a V?-” ‘complementary’ to V”". 
If Bg refers to Vj", the fundamental tensor of V?-™ is C3} = 63—BZ, | 
and the Eulerian curvature tensor of V"'~™ is (—3). Thus the con- 
ditions for a conjugate congruence of order p in V?~™ are 


[C; p]’ 05 = 0 (r = 1.,...,p; « = l.,...,”). 
Since O48 = Q;* Bi g+ B; *O! 9 
and Bi gt = 0, 


these conditions reduce (in number) to 


[C; p] Of (*=8 (fe = 1.,..., 9; 6 = 1,..., 0). 


The case of m = 1 is of interest. If V! is composed of the con- 
1 
gruence A“ along the curve C, the conditions [C; | for a conjugate 
congruence of order p in the complementary V”-! are 
g vi 3 1 
Dr 


ZrTre 


I > 
i.e. BA (* = 0 @ = I....,9), (16) 
: 3.2 


where A* is the first normal congruence, and & the first curvature, 
2 1 


of the congruence A*. 
If A“ is a parallel congruence, every congruence in V~" is conjugate 
of order n—2; otherwise, conditions (16) are equivalent to 
{*A, = © (r = 2,...,p+1); (17) 
r 


i.e. {* must be perpendicular to the first » normal congruences of 
the given congruence. 

The (unique) conjugate congruence of highest order is the con- 
gruence A“. 


If A* is taken to be tangential to C, then for (* to be conjugate | 
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of order p in the V?-1 normal to the curve, it is necessary and suffi- 
cient that (* be normal to the osculating (p+1)-space of the curve 
at each point. 

REFERENCES 


. E. Bortolotti, Rend. Ist. Lombardo, 64 (1931), 441-63. 
». L. P. Eisenhart, Riemannian Geometry. (Princeton, 1926.) 
. H. A. Hayden, Proc. London Math. Soc. (2) 33 (1931), 22-31. 


. , ibid. (2) 34 (1932), 27-50. 
. V. Hlavaty, Rend. Acc. Lincei, (6*) 12 (1930), 647-54. 
A. J. McConnell, Proc. London Math. Soc. (2) 28 (1929), 501-9. 
J. A. Schouten, Math. Zeitschrift, 30 (1929), 149-72. 
J. A. Schouten und E. R. van Kampen, Math. Annalen, 103 (1930), 752-83. 
. D. J. Struik, Grundziige der mehrdimensionalen Differentialgeometrie .. . 
(Berlin, 1922.) 














AN EXTENSION OF THE METHOD OF 
GENERALIZED PLANE STRESS 


By W. M. SHEPHERD (Bangor, N.Wales) 
[Received 12 December 1932] 


In this paper the method of generalized plane stress is extended to 
apply to plates of which the thickness decreases uniformly in one 
direction. 

1. The stress function 


! 


We use cylindrical coordinates r, 6, z and take the planes @ = +a 
to define the two sides of the plate. 

Following a method analogous to that used in generalized plane 
stress,* let us assume that 


~ 


6¢ = 0 throughout, 


bz = 6r = 0 when 0 = +a. 
Let 


x x Ta 
] ~ - a ~ 1 ~ 
rr = — | rr dé, zz = — | 2zd0 m = — | rzd@. (1.1) 
2a , 2a , 2a 


The body-stress equations, in the case of equilibrium with no body 
force, aret 


orr 1 érd. erz . rr—00 


eon he reer =e 

or |r 06 az ' r 

érd 1000, Oz | 2rd _ (1.2) 
Or  v OF a2 ¢ 


Orz 1002 . 022. rz 





Orr 00 oe fr 
On integrating the first and third of these with respect to @ from 
@ = —a to 0 = a, we get 


—(rrr) + —(rrz) 0 

cor C2 

A (1.3) 
— (rrz) + — (rz) = 0 

or oz 


* Love, Mathematical Theory of Elasticity, 207. (1927.) 
t Ibid., p. 90. 
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The second of equations (1.2) gives no information about rr, 22, or 72. 
The equations (1.3) are satisfied if the mean stresses are derived from 
a stress function x by the equations 
"x 2 = x, in o. (1.4) 
r Oréz 


a 1 
rr —s n.9? 
r 02? r Or 


2. The differential equation satisfied by the stress function 

To find the differential equation satisfied by the stress function, 
we use the stress-strain relations. 

The mean displacements are defined by the equations 


U, = = [ u,d0, U,=> [ u, dd. (2.1) 


r , 
2a , 2a J 


x x 
Then precisely as in the theory of generalized plane stress* for a plate 
of constant thickness, we obtain the equations 


where 1—o = 1/(1+7) and 7 is Poisson’s ratio. 


If we denote ox, ox by #y, the equations (2.21), (2.22), (2. 
é dz? 


9 


give, with (1.4), 


¢ ou, 

— = 
Cz 

(2 | au) _ 1 ay 


C2 or Tr Oorcez 


These equations must be consistent, and so, by eliminating U, and 
U_. we obtain the differential equation satisfied by x, which is found 


to be o ne 
92/1 92, _ 21+) ex (2.4) 


2 


lr “J rs Gz 


* Coker and Filon, Photo-Elasticity, 129. (Camb. 1932.) 
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3. The displacements 
32) we put 1 52, = 7. 
r dz? 
and integrate with respect to z, we get 
la 
U, = (1—o) 1 2X 4 Fy, (3.21) 
dz =r Oz 


where f(r) is an arbitrary function 
To obtain U,, substitute for U, from (3.21) in (2.33) and integrate 


We then obtain 
Li 
—(1—o) 1X _xX_ p(n) 
dr ror fr 


with respect to 
9 = 
2uU, = 

where g’(r) is an arbitrary function, written for convenience in this 


form. 
ry eo 
hence we have the further condition 


(1-0) 8% = =* —2f"(r)—-g"(r), 


are g'(r) is & 
These values of the displacements must also satisfy (2.31), and 


(1—o) 96+ (1+ ){zf(r)+-9(r)}] 


or 
b= $+(1+nXef(r)+9(r)} 


If we write 
we then have the displacements given by 


ul, 


where ¢ is defined by 
5 — 2(1+ mx. 3.5 
Y= r3 ( 5) 


and 
It is easily shown that a function 4% can be found which satisfies 
In this way % is defined except for terms of the 


(3.4) and (3.5). is way 
form ar+bz+crz, which correspond to states of zero mean stress 
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4. The solution of the differential equation for the 
stress function 
4.1. Solutions independent of the elastic constants. 


If x is independent of 7, we must have 


a =0 (4.11) 


and oe ail (4.12) 
satisfied simultaneously. 
Equation (4.12) gives x = R,+Rpz, (4.13) 
where R, and R, are functions of r only. 
On substituting in (4.11) we find that R, and R, satisfy the 
equation d2 (1d2R 
war 7A pits 
The solution of this is R = of*+br*+er, (4.14) 
where a, b, and ¢ are constants. No constant term is included as it 
gives rise to no mean stresses. 
We then obtain the solutions 
xX, = art+br', (4.15) 
no term in cr being included, as it gives rise to no mean stresses, and 
X2 = (ar*+6r3+cr)z. (4.16) 
These solutions are discussed in § 5. 
4.2. Solutions in the form of power series in r. 
Let us assume that Ee Mf 
x= > han (4.21) 
n=0 - 
where the Z,, are functions of z only. 
On substituting this value of y in (2.4) and equating the coefficients 
of powers of r to zero, we obtain the results 
Zo = az*?+62', 
Z, = c+dz, 
Z, = y(2a+6bz), 
20,.1(n+2)*—(1+n)], 20% _ 9 aaa) 


Anis . 
n+3 (n+-1)(n+-2)(n+3) n+ 
where n = 0, I, 2.,.... 
3695.4 
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It follows that Z, and Z, are arbitrary and independent of Z,, Z,, 
and Z,. The term (c+dz)r is included in the solution already obtained 
in (4.1) and so is omitted. 
We then get a solution 
xX = az*+bz3+ n(a+ 3bz)r?, 


or, since @ and 6b are independent, 


Xs = 2+ 77, (4.23) 
Xe = 2+-3y2r*. (4.24) 


The remaining terms give a solution 


A yr 

i> 4 

x= > Z,>> (4.25) 
n 3 , 


where Z, and Z, are arbitrary and the succeeding terms are given 
by (4.22). The question of the convergence of this series is not con- 
sidered, as only terminating series will be used in this paper. As 


examples, let us take some simple values for Z, and Z,. 


If Z, = 602", B, = O, 
then Xs = 10r°z?—(3—n)r’. 

If Z, = 602%, Z,= 90, 
then Xe = 10r?z8—3(3—n)r*z. 

If Z, = 0, Z, = 10802?, 
then X7 = 45r4z?—(8—n)r®. 

If Z,= 0, Z, = 3602', 
then Xs = 15rtz?—(8—n)r*®z. 


4.3. Solutions involving logarithms. 

We have obtained solutions involving the two arbitrary functions 
Z, and Z,, but since the equation is of the fourth order, we expect 
to find solutions involving four arbitrary functions. 

By analogy with the method of solution of ordinary differential 
equations in series, we assume 


0 
x —— Xo log r= Zz, a r, 
0 


n 
where x, is a solution of (2.4) obtained from the functions Z, and 
Z,, and the ¢,, are functions of z only. 
We then substitute in (2.4) and evaluate the ¢,, by equating the 


coefficients of powers of r to zero. 
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As an example of this type of solution we find that if 


2 i Xo = Xs 
another solution is 


67—2n—39n? 5 

3 T 28y3— — > 
2(3—m)(1+-) 2(1+-m) 
Further solutions, corresponding to yx., x;, and xg, may easily be 
obtained by this method. 


Xo = logrx»+ zr. 


5. Applications of the solutions 
5.1. Simple tension, Euler-Bernoulli flexure, De Saint Venant 
flexure. 
In equation (4.15) we have the solution 


xnu= art +-br°. (5.11) 
The stresses are 


~ 


iv = 9, r 


( 


a= @, zz = 12ar+-6b. (5.12) 


x 


The corresponding value of ¢% is 
wb = 3(2ar+b)z?—2ar?— 3br?+ 4(1+ n)(ar3+-3br?). 
[t follows that the mean displacements are given by 
Pia > A-1s 2] 4-2) 1 

U, = —6E-Ya(z?+ nr*)+ br} | (5.13) 
T > 7 ' ? ° 

U, = 6E-(2ar+b)z J 
where # is Young’s modulus. 

When applied to a strip with edges r = 8 and r = y this solution 
corresponds to uniform tension and to bending by terminal couples, 
and so is of no special interest. 

Now, in equation (4.16), we have the solution 


Xo = (ar*+bre+er)z. (5.14) 
The stresses are —— 
rz = —(4ar?+3br+c/r) }. (5.15) 
2z = (12ar+6b)z 


The function ¢ is 
yb = (2ar+b)z3+-2{(4ar3+br?— 2c log r)(1+ n)—2ar?— 3br*}, 
and hence the displacements are 
U, —E 1f2a23-4 6yrz(ar-+-b)} 
U, = E- 1f[ (2ar-+-b)32?— 2ar3— 3br?|— (1+ »)[Zar3+-c(1+-2 log r)]}}. 
(5.16) 
K2 
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This solution may be applied to a strip with edges r = 8 and r = y. 
The constants a, b, ¢ can be so chosen as to make these edges free 
from traction. We then have the solution for De Saint Venant flexure 
or bending by a transverse load at one end. 

The distribution of stress is different from that in a strip of uniform 
thickness owing to the term c/r in rz, but 2 is of the same form. 


5.2. A strip bounded by r = a and r = b and under normal pressure 
on these edges. 
In equation (4.23) we have the solution 


an ye KO 
Xs = 2°-+- 77". .21) 
The stresses are 


rT = 


If b = 0, we have the case of a knife pressed by a uniform load 
on the back r = a against a plane object to be cut. There are, how- 
ever, tractions across the ends which we may suppose to be z = -Le. 
It is not possible to annul these tractions completely by any of the 
stress functions found, but we can superpose a system derived from 
x, Which will make the unwanted tractions at each end form a system 
in statical equilibrium. The effect of such a system, except near the 
ends, may be neglected. 

The tractions over the end z = c, due to x5, are equivalent to a 
force 4aay in the line r = 0 and a couple of moment* 2aa?7. 

The stress function required to annul this force and couple is 

x = (nr*—2nar*)/a?. 
The complete stress function is 
xX = 22+ nr?+ (nrt/a*)—(2nr3/a), 


and the stresses are 


72 Q@, Zz 


> 


2 
, 


wage 2n 12yr_ 12y 
y 


a> a 


5.3. A strip bounded by r = 0, r = a under traction on r = a. 

The solutions of the type (4.25) are applicable to this problem since 
all give zero traction on r = 0. 

First consider a beam under uniform normal load on r = a and 
supported by shearing stress across the ends z = +e. 


* This is only approximate. An exact value could only be obtained if the 
complete, instead of mean, distribution of stress were known. 
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The stresses due to the stress function 

X = X1—6ax, (5.31) 
oe rr = —30(4ar2—3r°) 
zz = —30[12az?— 18rz?+-4(3— n)ar?+ (8—n)r>] (5.32) 
rz = 360zr(a—r) 

Onr=a rr = —30a?, rz = 0. (5.33) 
We have then a uniform normal load on r = a, and at z = +c we 
have a set of normal tractions and shears. The shears at z = ¢ are 
statically equivalent to a force of magnitude* 60aca* in the direction 
of r increasing. This force is equal and opposite to half the load on 
ry =a. The normal tractions may, as in (5.23), be combined with 
a set of tractions derived from x, to give a system in statical equi- 


librium. 
We can of course combine the solutions of (5.2) and (5.3) to give 
that for a knife pressed by transverse forces at the ends against the 


object to be cut. 


5.4. A uniformly loaded cantilever. 

The solution of (5.3) can be applied to the case of a cantilever free 
at the end z = 0 and clamped at the end z = c, provided that we 
again add suitable terms to the stress function to annul the force 


and couple resultants across z = 0. 


5.5. Other soluble problems. 

3y suitably choosing Z, and Z,, we can obtain the stress functions 
corresponding to other distributions of tractions on the edge r = a 
when the edge r = 0 is free from traction. 

If we require either a distribution of normal traction on the edge 
r = 0 proportional to z* or a higher power, or of tangential traction 
proportional to z or a higher power, it is necessary to use the stress 
functions involving logr, of which x, is an example. 

In order to find the complete solution for any system of traction 
on the two edges of such a strip, it would be necessary to find the 
solutions of (2.4) in the form F(r)cos mz. The functions F(r) so found 
are, however, so unmanageable as to make such an investigation 
unprofitable. 


* This is again only approximate. 





ON THE DISCONTINUOUS INTEGRAL OF 
WEBER AND SCHAFHEITLIN 
By E. T. COPSON (St. Andrews) 
[Received 15 February 1933] 
1. Introduction 
THe Weber-Schafheitlin integral, namely, 


Iz { J, (t)I,(utyt> dt, (1.1) 


where the parameter wu is positive, has been the subject of many 
papers.* Quite recently, Dixon and Ferrart have shown how the 
discontinuity of this integral, according as u >, =, or < 1, may be 
clearly demonstrated by applying Cauchy’s theorem on contour 
integration to a similar integral involving Bessel functions of imagi- 
nary argument. 

In the present paper, a new method of evaluating the integral (1.1) 
is given. From this, it appears that its discontinuity, whether in 
actual value or in the form of expressing that value, as the case may 
be, is directly dependent on the discontinuity in value of 


wi 


. 
| efxy—% dx, 


a 
according as @ >, =, or < 0. In fact, when there is discontinuity 
in the actual value of J, namely, when —1 < R(A) < 0 and p—v is 
an odd integer, the two discontinuities are of precisely the same 
nature. 
2. The fundamental formulae 
When uw is positive, the integral J in (1.1) converges provided that 
R(iut+v+1) > RA) > —1 (u # 1), 
R(u+v+1) > R(A) > 0 (wu == 2): 
We shall find it necessary to restrict still further the parameters 
A, », v in order that certain transformations of integrals may be 
carried out, but shall appeal to the theory of analytical continuation 
to show that the values of J so obtained are valid whenever J con- 
verges. 
* See, for example, Watson, Bessel Functions, 398-410. (Cambridge 1922.) 
+ Quart. J. of Math. (Oxford), 1 (1930), 122—45. 
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We replace J,(ut) in I by a modified form of Hankel’s contour 
integral, namely, en 
1 ¢t\ t?\ dx 
J (ut) = exp/ «*z — — }|—__, 
r(ul) 2art | (=) r{ z\en 


cC—RI 
where the path of integration is the straight line R(x) = c > 0; this 
formula is valid when R(v+1) > 0. Next, we invert the order of 
integration, to obtain 


eB ewe , 
[ dx picoceey | v-rJ, (tes dt. 

c—wi 0 

Sufficient conditions for the validity of this process are that 
R(wt+v+1) > RA) > F, Riv) > 0, 
which ensure that the repeated integral converges absolutely. 
But, under these same conditions, the integral with respect to ¢ 

may be evaluated by means of Hankel’s generalization* of Weber’s 
first exponential integral, so that 


Vi a ie | | i de we oe sr, 
[== ete | —_— iF, — A atl;—2 dz 
2u’T(u+1) J aheti+A-p) 2 


c—@® 


(2.1) 


C+ 
e(u—lx EF A+p+1—v 
giettA—p)r 2 


Cc—Ri 


sm+sa) dx, 


(2.2) 
the second of these integrals being derived from the first by Kummer’s 
transformation.t| The formulae (2.1), (2.2) are fundamental in our 
subsequent work. 

Now, by Barnes’s asymptotic expansiont{ of the generalized hyper- 
geometric function, we see that, when 2 = c-+-iy, where y is large 
and positive, 

\(du+$v+}—DAje* 1—A 
Putt h— ews p (utytI—A, 
Aw’ DT (w+ Late r+t—-w 2 
ae i(v—A—p—Do(u? x Ti 1 i Loe ewan 
oe ) 5 =) . (s¢+3v+9 24) ; . (2.3) 
wralt I\($A+ 34+ 43—}v)24w’al +” 
hence the integral (2.1) converges when R(u+v+1) > R(A) > —1, 
R(v) > —1, (wu £1) or when R(u+v+1) > RIA) > 0, RW) > —1, 
* See Watson, loc. cit., 393, (2). 
+ Viz. ,Fy(a; p; x) = e*,F\(p—a;p; —x). See Watson, loc. cit., 102. 
t Cambridge Phil. Trans. 20 (1908), 259. 








136 E. T. COPSON 
(wu = 1).’ It can be shown, by analytical continuation, that the equa- 
tions (2.1) and (2.2) still hold under these less restrictive conditions. 

The asymptotic formula (2.3) also shows that, when wu = 1, 
R(A) = 0, the integral in (2.1) oscillates finitely, save when A = 0, 
though in the latter case it exists as a Cauchy principal value. 
3. The general case 

In certain particular cases, the generalized hypergeometric func- 
tion which occurs in (2.2) is a polynomial, a fact which then simplifies 
the evaluation of the integral. In the general case, we use the 
formula 

1( as lpr x) LF x; p32) ex(1—¢p- x—1(] —t)*- dt 
I'(p) ; 

which holds when this integral converges. 

If we suppose, temporarily,* that 

R(u+v+1) > RA) > | R(u—v+1)], 
we obtain in this way from (2.2) the equation 
1 sa 
iw’ T(SA+4u-+4—4y) © 
ctat 1 


a a{u?—t) 
x ae | : thu +v-A-1)(] _¢)}A+u—v—)) dt 


I= 


ah tA+1—p) 
c—ai 0 


1 
= DapiwT(GA+du+3—h) * 
1 c+ wi 


F (ut 
x | a te+v—A-1)(] - {)8A+u v—1) ies Hips 


dx, (3.1) 


giv +A+1—p) 


0 c—ai 


the inversion of order being justified by the absolute convergence of 
the repeated integral. 
Now, since R(v-+A—p—1) > 0, we havet 


C+ 
ex(u*—t) (2 2__g)Kv+A p—1) 
dx = aan 


or 0, 
3H) 


av +A+1—p) 
cat 


according as u2 > or < #¢; it is from this fact, as we shall now show, 
that the discontinuity in form of the value of J arises. 
* We could have avoided these restrictions by the use of the Pochhammer 


double-circuit integral for ,F,. 
+ Cf. Whittaker and Watson, Modern Analysis, 246, Ex. 1. (Cambridge 


1920.) 








ON THE INTEGRAL OF WEBER AND SCHAFHEITLIN 137 
When u > 1, the inequality u* > ¢ is satisfied over the whole range 
of integration with respect to ¢, so that, in this case, 
1 
—————- x 
j— Prt PtI— Ip) 


x [ (tt) nngte tay aise (3.2) 
0 
: ur-B- 1D (4u+4v+43—}A) 
AD + A+4— $y) P41) 
. p(ttet io ee thigh (3.3) 


9 ’ 9 


by a well-known integral* representing the hypergeometric function. 
The restrictions under which (3.3) has been proved may now be 
relaxed by the use of the theory of analytical continuation; in fact 
(3.3) holds for w > 1, whenever J converges. 

On the other hand, when u < 1, the inequality u* > t is not satis- 
fied over the whole range of integration with respect to t, so that 


we now have 


XxX 


] 
2 wD (3+ $u+-4— dw) (d+ A+ $— $y) 


—t) +A—p—Dgh(u +v—A-1)(] —4)8A+p—v-1) dt. (3.4) 


When uw < 1, we make the transformation ¢ = ru? in (3.4) and obtain 
wh(et++I—-P) 
AT (A+ d+ 4— WT O+1) 
tedtk ¢i-1—{— 
x F(t! | v+l1- Av l = Ht lyw', (3.5) 


9 


a result which still subsists whenever J converges. 

Finally, when w= 1, the integral in (3.4) is a case of the first 
Eulerian integral, so that 

PAM Qu tiw+si—P) 3.6) 
AD ($A+ du+4— $v) P(4v+ A+ $—fu) (A+ Sut vt hy)’ 

In this case, too, we may show, by analytical continuation, that 
(3.6) still represents J, when the only restrictions on the parameters 
are that J should converge. 


* Whittaker and Watson, loc. cit., 293, Ex. 1. 
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Now when R(u+v+1) > R(A) > 0, J converges for all positive 
values of u, and a comparison of (3.3), (3.5), (3.6) shows that it is 
continuous at uw — 1. The form of the value of J is, however, discon- 
tinuous, since the value for u > 1 is not the analytical continuation 
of the value when u < 1. There is, however, a genuine discontinuity 
when —1 < R(A) < 0, since the integral does not then converge for 
all values of wu. 

4. The critical case 

It follows from the asymptotic expansion of the Bessel function 
that the integral (1.1) converges for all positive values of wu, when 
v—p is an odd integer, provided that A lies in the strip 

—1< R(A) < R(wt+v+1); 
further, for any fixed value of w, the integral is an analytic function 
of A, regular in this strip. 

Now in this particular case, the expression (3.6) is also a regular 
function of A in the same strip, the singularity at A = 0 being a 
removable one. Thus, by analytical continuation, the formula (3.6), 
or its limiting form as A — 0, gives the value of J, when v—p is an 
odd integer, in the extended range —1 < R(A) < 0. 

But when we regard the Weber-Schafheitlin integral as a function 
of u, the discontinuity of form becomes an actual discontinuity in 
this critical case. Although this discontinuity may be demonstrated* 
on the lines indicated above, it seems to be of some interest to show 
how it arises by the direct method of the present paper. 

Let us consider, then, the integral 


J= [v (ut) dt, (4.1) 


a-p-i 


(t)J,, +p 
0 

where p is a positive integer or zero; it converges when R(a) > 0. 

It follows from (2.2) that 





c+ot 
a? (u*—1)x 
J = —. (a) - S \(—p;a—p;2) dx 

2ri u**PT'(a—p) git 

c—ni 
C+at > 
— I(x). p! we [ e(u?—1x : I\(l—a+r) dx, 
2riu*+?T\(a—p)(1—a+p) r!(p—r)! al 


c— wi ia (4.2) 
the integral being a Cauchy principal value when u = 1. 


* See Dixon and Ferrar, loc. cit., 132. 
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Now, when u? < 1 and s is a positive integer, 


c+ai 
(u?—1)x 
[ : —dz = 0, 
a 


J 
C— at 


so that J vanishes when u < 1. But when u = 1, 
c+Ri 
J = ES . D(a)P(—a) lim | dx _ 4(—1). 

2nt T'\(a—p)T(1—a+p) r+» x “ 
c-Ri 
> 1, we have,* by (2.1), 
TU (o+ ) 2 

(«) J : F(a; a—p; —x) dx 


2riutt?T\(a—p) J ai+P 


Finally, when u? > 
F as 


(—2x) 


ees | “es T(a+r) o* 
eS oe 4 T(a—p+r) r! 


2m ust? 


F(a) ] 
= F\ a, —p;«a—p;—]. 
u*-?T'(a—p).p! (« —s ) 
We may, then, sum up the results of this section in the formula 
‘ dn Fla, —p3a— Pia 4 , 
ur? T(a—p). p! 
| Jo rp(Ut)I s(t) o = 4(—1)? 


0 
0, 


or < 1, provided that R(a) > 0. 


according as u >, = 
* It should be observed that this change of contour may be applied to (2.1) 


in the general case and that it leads directly to (3.3). 

















DU BOIS REYMOND’S CONSTANTS 
By G. N. WATSON (Birmingham) 
[Received 6 April 1933] 
THE constants which I shall discuss in this paper have their origin 
in the following theorem: 
Let s a,, be an oscillating series whose upper and lower limits, U and 
1 


n 


L, are finite. Let 4(t) be a function with a differential coefficient which 
is finite when t > 0, such that 





, * |dd(t)| 
lim ¢(t) = 1, | |‘ H(t) dt < o. 
t>+0 J | dt 
2 0 
Then, if the series > a, (nt) converges whenever t > 0 and has the sum 
n=1 


S(t), the upper and lower limits of S(t) as t > +-0 lie between the two 
numbers = 
\(U+L)44U—L) | 


0 





a4(t)| 
di | dt. 


The history of the theorem is as follows: in Du Bois Reymond’s 
work* on trigonometrical series, the series > a,, was a trigonometrical 
series with a, = 0(1) and ¢(¢) had the special value ¢-*sin*t. Du Bois 
Reymond obtained the upper and lower limits 


}(U+L)+4(U—L)(1+-), 


9 9 
where e< 2+-+-, 
7 a 


by an ad hoc argument. The type of argument which he used was 


apparently sufficient? to prove that 


and the latter constant had, in fact, appeared in 1854 in Riemann’s 
Habilitationsschrift,{ where it was assumed that > a, was convergent 
so that U = L. 

* Abh. der Kgl. Bayerischen Akad. zu Miinchen, 12 (1876), 119-66. 


12 
+ Hobson, Functions of a Real Variable (1907), 732-4. 
t See his Ges. Math. Werke (1892), 227-64 (247). 
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Later, by an argument which is applicable to a general function 
¢(t), G. C. Young* proved that, when ¢4(t) = ¢-*sin*#, 


is 2) 
d (sin*t ee 
l+c = = , dt <1+-—-, 
J jdt\ # 7 
( 
and the theorem in its general form was actually given by Hobson.f 
Since Du Bois Reymond’s work exhibits the importance of the 
numerical value of c, it seems appropriate to associate his name, 
rather than Riemann’s,{ with such a constant, and to describe 


2) 


[ > |*— 
dt 


0 
as Du Bois Reymond’s constant associated with ¢(¢). 
In recent work on trigonometrical series,§ derived series in which 


: sin t\* sin t\4 
#(t) has been given the values te and 5 have been the 


subject of investigation. It consequently seems worth while to 
examine the constants c,, defined by the formula 


wn 


ca = fie (em9| a 
» Jee 


) 
[ shall show that c,,, is expressible as a polynomial of degree m in 
e? with rational coefficients.|| There seems to be no similar simple 
formula for ¢y,,,,, and the sum total of the results of my analytical 
attacks on cz has yielded merely an integral which does not lend 
itself readily to numerical computation plus a series whose sum is 
fairly easy to compute; I have, however, found it easy to compute 


c, by elementary methods. 
sin f\2" 
i 6 = ("PY 
t 
it is easy to see that the minima of ¢(¢) in the interval (0,00) occur 


* Messenger of Math. 49 (1920), 73-8. 
+ Functions of a Real Variable, 2 (1926), 221-4. 
{ Obviously all that matters in Riemann’s work is that c should be /inite; 
its numerical value is irrelevant. 
§ Verblunsky, Proc. London Math. Soc. (2) 31 (1930), 387-406; Journal 
London Math. Soc. 6 (1931), 82-4. 
The corresponding coefficients in the expression for 2°"~'(m—1)!c,,, are 


integers. 
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at t = 7, 27, 3z,... and the maxima occur at t = 0, w,, we, @3,..., 
where w,, @2, @3,... are the positive roots of the equation 


sint —tcost = 0. 


Also b(n) = 0 in =. 25), 
] 
¢(a,,) = cos*"a,, = (a = 1,3,3.,..2. 
(1+ o>)” 
Hence 


1+-Com _ {6(0)—¢(z)}+- |(7)— $(a) a 
+{$(@)—$(27)}+ |b(27)—$(w,)|+-..- 


“ l 
so that Con, = 2 MS al ym 
+ Wn)" 


n=1 


Now consider the integral 


1 zsin z dz 
2ni J sinz —zcosz (1+22)” 


where the path of integration is a (large) circle with centre at the 
origin. As the radius of the circle tends (discontinuously) to infinity 
in such a way as to avoid the poles of the integrand, the integral 
tends to zero. 

The poles of the integrand are at the origin, where the residue is 3, 
at the points +a, (n = 1, 2,3,...), where the residues are (1+ w7)~”, 
and at the points +7. 

By putting z = +7¢, we see that the residues at +7 are both equal 
to the residue of 

Csinh f l 
~ sinh £ —fcosh ¢ (1—@2)™ 


at € = 1, that is, they are equal to 
(—)* [ees Csinh a. Oe 
(m—1)!| dg"—|C cosh £ — sinh ¢ (C+-1)") Jr.y 
Hence, by Cauchy’s theorem, 


— (—)"2 [ d- / {sinh € > 
(m—1)!| df"-1|f cosh £ — sinh Z (€-+-1)”J fa. 


—— 
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Csinh f hy f(et—e-) 
Ccosh£—sinh£ (£+1)e-$+-(f— 


= Uet—et) » Shee 


1)eé 


and therefore 
Csinh f 1 . 
fcosh{—sinhZ(f+1)™ —— (€+1)"™42 ((41)"' 
9 l 


__\r-1 _y)r-1 ere] l ea = \. 
PE "a Cet eel 


Hence 
(a ii | Csinh Z 
df" |C cosh f — sinh ¢ (Z 
eee 


_ (2m—1)! 


2297)! 2m—1(m—I)! - 
é 1 
¢=1 


. qm —r gort| 1 
= (f+ 1)er+t (¢+1) ym" ir | (¢+1) m+r— 1 


dtm —r 
m—r 


(2m—2)! 4 > (m—1)! | >, (Ym  C,(2r)8e" x 


m—r)! 


22mm! 

s=0 
/{ (2m—s)! (2m—s—1)! | (2m—s—2)! 
~~ 92m—8(m+r—1)! * 22"-8-(m+r—2)!J | 


* | 9% 2m—8+1(m +r)! 


That is to say, 
3 (2m—2)! 

~ — 22m—-1m!(m—1)! ' 

~~ ~< — (—)*(4r)8e?" {(2m—s)! (2m—s—1)! , ,(2m—s—2)!) 

Ly Lx 22m 1s!(m—r—s)!| (m-+r)! (m+r—1)! (m+r—2)! J’ 

rT s=( 


which is the required expression for C.,,. 


This formula gives the following results: 
P Co = 0-19452 80495, 


2c, = e*—7, 
. = 0-00524 07047, 


SC, 2 e4— 4e2— 25, 
32c, = e&—6e4+ 3e?— 98, C, = 0-00022 06747. 
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The value of c, should be compared with the value of G. C. Young’s 
approximation 27-7 = 0-20264. 

Up to a point the investigation of c, is very similar to that of ¢,,,,. 
We observe that the maxima of ¢-*sin*¢ in (0,00) are at t =0, wp, 


@,,..., and the minima are at w,, @s3,..., Whence 


5‘ sina, or l 
>, (Ya = 2 > aa 
7, (1+ o;,)? 


n=1 os n=1 


Now consider the integral 


ee i sin®z dz 
27i | (sinz — zcosz)cosz (a?+-2?) ee 
where the path of integration consists of two (large) equal semi- 
circles on the left and right of the imaginary axis respectively, both 
having their centres at the origin, the ends of the semicircles being 
joined by loops to exclude the points -+-ai from the interior of the 
contour. 

The poles of the integrand are at the origin, where the residue is 
3/a, at the points +7a,, where the residues are (a?+@;7)-?, and the 
points +-(n—4)z, where the residues are —{a?+-(n—}4)?x?}-4. 

It follows from Cauchy’s theorem that 


319 ~ ] ] aa ] | 
at? 2 leew) En DP 


1 n 


(ai—) (—ai—) ? 
sin?z dz 


(sin z —z Cos z)COs z (a?+2*) 


We sinh?¢ dt 
} J (€cosh f — sinh £)cosh ¢ ,/(a?—Z?) 


sinh2x dx 
_ » 9 9 
(a cosh 2 — sinh x)cosh x 4/(x?—a?*) 
4 sinh?ax dx 


2 
T | (ax cosh ax — sinh ax)cosh ax ,/(«*—1) 
1 


Now differentiate with respect to a, and put a = | after differentia- 
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tion; we get 
3+2 S . —§ . 
(1+)! + {1+ (n—4)?x}} 


i a) 


[ d sinh?ax dx 
a J | da(axcoshax — sinhax)cosh ax },_, ./(z?—1) 
1 


2 [ d { tanhx x dx 
7 J dx\x—tanha} ./(z?—1) 
1 


[ have not discovered any simple method of computing the integral 
on the right. 


It is possible to compute c, directly in the following manner. We 
have* mw, = 4-49340 94579, @, = 7°72525 18369, 


while, if (n+ 3)7 = a,, we have, by a formula due to Euler, 
1 2 13 146 
3a3 1505 10507” 


l . 19 323 19519 
whence — —- a. f 


— - —-+- — - — 
2\% ‘ o> * Gat * Oe * ennG. itt * 
1+)! o8 © 205 © 807 800%  2688a!; 
From this result it is easy to compute 


2 


1 
>. Fey 


n=3 “s 
with the help of Glaisher’s tables of 


S l 
Qn 1. 1\m" 
. (2n-+1) 


This procedure gives ee 
8 procedure give Cz = 0-02825 17642. 


Finally I remark that 


we . 3 
- + sin’a, 5 
an sso. 
n=1 


* The values of w, and aw, to four places of decimals are quoted from 
Lommel by Jahnke u. Emde, Funktionentafeln (1923), 3. I computed the values 
stated directly by interpolation in Andoyer’s fifteen-figure trigonometrical 
tables. 
+ Introductio in analysin infinitorum, 2 (1748), 319; this is a special case 
of Stokes’s formula for zeros of Bessel functions applied to J;(2). 
Quart. J. of Math. 45 (1914), 141-58. 


3695.4 


L 
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is easily evaluated by considering 
1 f¢ 2? dz 


2mt J sinz —zcosz (1+2?) 


9? 


and it is found to be equal to 
3—4(e3—3e) = 0-01732 71405. 

While computing c,, I also computed this constant by the same 
method and obtained a value differing from the correct value by 
a unit in the eleventh place of decimals. 

It is, of course, possible to compute c;, ¢;,... in the way in which 
[ computed c,, and it is also possible to obtain integral representa- 
tions of c;, c,,... analogous to the integral for c,; given above. 

Also, by considering the contour integral 
:.F 2* dz 


2771 . 





sinz—zcosz (1-+-22)mt!’ 


it is possible to prove that 


1 x 


wD . 9 
) 3 gin?*+1 gy eq (—)""! 
* -2m+1 Bi } 2\m+4 
en ( l +am7,) she 
n=1 n=1 


is expressible as a polynomial of degree 2m+1 in e with rational 


coefficients. 




















ON CRITERIA FOR UNIVERSALITY OF 
TERNARY QUADRATIC FORMS 
By ARNOLD E. ROSS (Pasadena) 
[Received 8 December 1932] 

1. In this paper we derive necessary and sufficient conditions that 
an indefinite ternary quadratic form 

Y = ao2?+ By?+ yz*+ pyz+oxz+rary (1.1) 
with integral coefficients «,..., 7 should be universal, i.e. should repre- 
sent all positive and all negative integers. We shall also show very 
simply that every universal form represents zero properly. If 
p, o, tT are all even, form V is called classic ; if at least one of p, o, 7 
is odd, it is called non-classic. If (1.1) is a properly primitive classic 
form, conditions for its universality are expressed in terms of its 
generic characters. If (1.1) is a non-classic form these conditions are 
given in terms of the generic characters of an improperly. primitive 
classic form f = 2. 

Universal ternary quadratic forms have been studied by Dickson* 
and Oppenheim.+ The present writer noticed{ that if the deter- 
minant of a universal classic ternary form is odd or a double of an 
odd integer, and its first invariant Q is equal to +1, then Meyer’s§ 
criterion for equivalence yields at once necessary and sufficient con- 
ditions for universality of a classic form VY in terms of its generic 
characters. 

[In the present paper we study both classic and non-classic forms 
and obtain the desired necessary and sufficient conditions without 


the use of the above-mentioned deep and difficult theorem of Meyer. || 
THEOREM 1. Let f be a properly primitive, indefinite, classic, ternary, 
quadratic form with reciprocal F and determinant D. The necessary 


and sufficient conditions that f be universal are 


D is odd or double an odd integer, Q = +1, (1.2) 
and, for every odd prime p dividing A = D, 
(F\p) = (—Q|p). (1.3) 
Dickson (2), (4), 21. + Oppenheim (6). 
2oss (7). § Meyer (5). 


A brief report on the results of this paper appeared in the Proceedings 
of the National Academy of Sciences, September 1932. There is a misprint in 
the statement of Theorem 2. 


L2 
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THEOREM 2. Let f, be a primitive, indefinite, non-classic, ternary, 
quadratic form. Consider an improperly primitive form f = 2f,. Let 
F be the reciprocal and Q, A the invariants of f. Then f, is universal 


if and only if f = 2f, satisfies the following conditions: 
O = +1, (1.4) 
the characters of f are 
(—1)?) = (—1)40-, (Fp) = (—Qip) (1.8) 
for every odd prime p dividing A = D, and, in case 4 divides A, 
(—1)i"-) = 1. (1.6) 


2. First lemma. In this section we reduce* the problem of 
representation of an integer m by a classic form (1.1), of determinant 
D with « equal to an odd prime or double such a prime, to the 
problem of representation of am by a form of type 

® = x?+ By, z), (2.1) 
where B(y,z) is a binary quadratic form of determinant aD. We 
show that if the form (2.1) represents a multiple am of m, then the 
original form ‘’ represents m. 


LemMMA 1. Consider a classic ternary quadratic form 


f = ax?+by?+cz?+ 2ryz+ 2sxz+ 2tay (2.2) 
with the invariants Q and A. Write 
ab—t? = OC, ac—s*? = OB, ar—st = —QR. (2.3) 


Multiplying both members of (2.2) by a, we get 

af = (ax+ty+sz)*+ (ab—l*)y?+ 2(ra—st)yz-+ (ac—s?)z*, (2.4) 
or, in notations (2.3), 

af = (ax+ty+sz)?+OQCy?— 20 Ryz+OQB2*. (2.5) 

Next write O(X,y,2z) = X?+0Cy?— 20 Ryz+OB2?. (2.6) 
Let a be a prime or a double of a prime not dividing determinant 
D=Q?A of f. Then (i) of f represents m, ® represents am, and 
(ii), conversely, of D represents am, then f represents m. 

Part (i) of the lemma is obvious. To prove part (ii) we assume 
that 4 
O(X,y,z) = am. (2.7) 

* Lemma | is a particular case of a more general theorem for quadratic 
forms in 7” variables and also of another theorem which reduces the problem 
to that of representation by a sum of squares with integral coefficients. Seo 
the above-mentioned report in the Proceedings. 
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By (2.6) and (2.3) 

O(X,y,z) = X?—fy?— 2styz—s*z? (mod a) 

= X?—(ty+<sz)? (moda). 
Since by (2.7) ® = 0 (moda), and since a@ is a prime or a double of 
a prime, either X —ty—sz = 0 or —X—ty—sz = 0 (moda). Hence 
there exists an integer x such that eX = ax+ty+<sz, where e = +1 
or —1. Substitution of this value of X into (2.7) and relation (2.5) 
yield the equality 
am = (ax+ty+sz)?+OCy?— 20 Ryz+OQO Bz? = af(x, y,z). 

Therefore m = f(x,y,z), and, since x, y, z are all integers, f repre- 
sents m. 

3. Second lemma. The proof of our second lemma employs an 
argument from the concluding part of Oppenheim’s* proof of his 
first theorem. 

LEMMA 2. Consider a form 

h = 274+ Ly? +2Myz+N2? = 2?+4(y,z) (3.11) 


of determinant 


d = LN—M?, where (L,2M,N) = 1. (3.12) 

If (L,d) = 1 and —L is a quadratic residue of d, (3.13) 

then h is equivalent to 2ary+ax?—dz?. (3.14) 
Multiplying both members of (3.11) by L, we get 

AL = La*?+v?+dz2?, v= Ly+ Mz. (3.21) 


By (3.121) and (3.131), —d is a quadratic residue of L. Also, by 
assumption (3.13), —Z is a quadratic residue of d. Thereforet there 


exist integers A, u, 7 such that 
LA? +-p?+-dr? = 0, My p (mod L) (3.22) 

with LA, », dy relatively prime in pairs. By (3.222) there exists an 
integer 5 such that 

p= 18+Mn,  (8,n) = 1. (3.23) 
Substitution of this value of » into (3.221) yields 

(8, 7) = —A?’, (6, ») = 1. (3.24) 

* Oppenheim (6), 182. 


+ By a theorem due to Legendre. For exposition see Dickson (3), Ch. VIIT 
(Statement on p. 130). See also Oppenheim (6), 179. 
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Hence there exists a transformation carrying ¢ into an equivalent 
form 

d(y,z) = (—A*, M,, ,), M?4+-N,»+* = M?—LN = —d. (3.25) 
Therefore f is equivalent to 

x2—)*y?+ 2M, yz+N, 2’, 
which by substitution of a for x—Ay is equivalent to 

x(a+ 2dy)+ 2M, yz+QN, 2”. 
Since (A, d) 1, it follows from (3.252) that (M,,A) = 1, so that 
there exist integers k, 1 such that Al—M,k = 1. Hence the trans- 
formation x, = Ax+M,z, 2, = kr+le 
is of determinant unity and yields the equivalent form 

22, y+ Hai+2F2,2,+ Gz. 
Dropping subscripts and replacing y by y+-ga— Fz, we obtain 

h ~ 2ay+ (H+ 29)x?—Hz?, 
where H = d since the determinant is unchanged. Hence, if £ is odd, 

h ~ 2xy+a?—dz? (3.26) 
by choice of g. If d is even, E is odd since h is properly primitive. 
Hence, to complete the proof of the lemma, it suffices to prove (3.26) 
for d odd and E even. In that case we choose g = —}E and see that 
h ~ 2xy—dz*. 

We replace « by «+4(1—d)y and get 2xy+-(1—d)y?—dz*. Next we 
replace « by «+dy+dz and z by y+z and get (3.26). 

4. Conditions (1.2) and (1.4).* If form (1.1) is primitive and 
non-classic, it represents all integers if and only if the improperly 
primitive classic form f = 2‘ represents all even integers. Therefore, 
henceforth we shall study a properly or improperly primitive classic 
form (2.2) and obtain necessary and sufficient conditions that it 
should represent all integers if properly primitive and all even 
integers if improperly primitive. Let Q, A be the invariants of f, 
and let F = Ax?+ By?+ Cz*+ 2 Ryz+28xz+2T ary (4.11) 
be its reciprocal. Multiplying both members of (2.2) by @ and C in 


succession and completing squares at each step, we obtain first (2.4), 


* See Oppenheim (6). 
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which becomes (2.5) in notations (2.3), and finally 
aCf = CX#?+QY?+O0Aaz?, (4.12) 
where X = ax+ty+sz, Y = Cy— Rz. (4.13) 
Next, in view of a well-known property of determinants, 
Aa+Tt+Ss=QA, Ta+Bt+Rs=0, Sa+Rt+Cs=0. (4.14) 
Without loss of generality we may assume that either a or, in case 
f is improperly primitive, $a is prime to 2Q0A. Then the coefficients 
C, R, B of the binary form 
(C,—R, B) = CY*?—2RYZ-+ B2?, BC— R? = Aa (4.15) 
have no odd prime factor in common. For if p divides all of C, R, B, 
it also divides both Aa, by (4.152), and the determinant A?Q of F. 
Hence p is prime to a and by (4.14) it divides S, 7’, A, which is 
impossible since F is primitive. If a is odd, we see in the same way 
that 2 also does not divide all of C, R, B, and hence (C,—R, B) is 
primitive. 
Let m be represented by f. Then (2.5) holds with f replaced by m. 
Hence for every odd prime p dividing Q, 


am = X* (mod p). 
If (m,p) = 1, then (m|p) = (a\p), (4.21) 


and therefore f does not represent }(p— 1) classes modulo p of integers 
m which do not satisfy (4.21). Hence if an odd prime p divides Q, 
neither f nor, in case f is improperly primitive, $f is universal. If f is 
improperly primitive, Q is odd and hence Q = +1. 

Next let f be properly primitive. Then a is odd and (C, — R, B) is 
primitive. If Q is divisible by 4, then, by (2.5), every odd integer 
m represented by f is congruent to a modulo 4, and hence f is not 
universal. Similarly is ruled out the case Q = +2, (C,—R, B) 
improperly primitive. Finally, consider a form f with Q = +2, 
(C, —R, B) properly primitive. Hence we may assume that C is odd. 
The equation (4.12) becomes 

aCf = CX*+ 2eY?+ 2eAaz?, 
where e = Q/|Q| = +1. 

If f, and hence aCf, represents all odd residues modulo 8, then 
y2+Aaz? = te(aCf—CX*), with X odd, represents all residues 
modulo 4. This latter is impossible if A is odd or is a naultiple of 4. 
It remains to consider A = 2A,, A, odd. 
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If f is even, so is X. Write X = 2X,. Then, by (4.22), if f repre- 
sents all even residues modulo 16, 


taCf = 2CX7+eY?+e2A, az 


represents all and hence all odd residues modulo 8. If $aCf is odd, 
so is Y, and CX?+ eA, az? represents all residues modulo 4, which is 
impossible since C and eA,a are both odd. 

Next assume that f is properly primitive and that its determinant 
Dis divisible by 4. By the aboveQ = +1, D = Q?A = A, and hence 
4 divides A. Then since (C,—R, B) is primitive, it is properly 
primitive by (4.152). Hence we may take C odd. By (4.12), 


aCf = CX?+QY? (mod 4), 


where CX?+-QY? has but three incongruent residues modulo 4 since 
both C and Q are odd. Therefore f is not universal. This proves that 
if f is universal, D is odd or double an odd integer. 

Summing up: If f is an improperly primitive form representing all 
even integers, its invariant Q = +1. If f is a properly primitive form 
representing all integers, then its invariant Q = +1 and its determinant 
D = QA = A is not divisible by 4. Since Q is odd, the reciprocal F 
of f is properly primitive. 


5. Choice of a and C. Up to the present we have been avoiding 
the use of the following normalization. 

Lemma 3.* Jf f is (i) a properly or (ii) an improperly primitive 
ternary quadratic form with a properly primitive reciprocal, then it is 
equivalent to a form f, such that (i) a, and C, are relatively prime and 
are prime to 2QA, or (ii) a,= 2x and « and C, are relatively prime 
and are prime to 2QA. Here a, is the leading coefficient of f,, and C, is 
the third coefficient of the reciprocal F, of fy. 

We shall outline a proof of it somewhat different from those 
referred to. 

Since F is properly primitive it represents properly an integer 
prime to any assigned integer and hence to 2QA. If QA is odd, then 
F represents properly an integer = 1 (mod 4) and one = 3 (mod 4). 
Hence we may assume that C is prime to 2QA and, if QA is odd, 


C = Q (mod 4). (5.11) 


* Smith (8), §§ 5 and 9. Dickson (4), 15-17. Bachman (1), 64. 
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Then the binary form 
yb = ax*+ 2try+by? (ab—? = QC) 
is primitive. 
Similar to (4.14) we have 
aS+tR+sC = 0, tS+bR+1rC =<, sS+rR+cC = QA, (5.12) 


The g.c.d. g of a, t, b divides ab—t? = QC and also the determinant 
D = Q°A of f. Since C is prime to 2QA, g divides Q and is prime to 
C. Then by (5.121) and (5.122) g divides sC and rC, and hence s 
and r. Therefore, since g divides Q, it divides cC and hence c, by 
(5.123). Thus g divides the g.c.d. of the coefficients a,..., t of f and, 
since f is primitive, g = 1. Hence ¢ is primitive. 

If f is properly primitive, so is 4. The determinant of f is 


abc-+- 2rst—ar?— bs?— ct? = Q?A. 


Let Q?A be even. If % were improperly primitive, then c would be 
even by the above and f would be improperly primitive. If QA is 
odd, QC = 1 (mod 4) by (5.11), whereas if % were improperly primi- 
tive, we would have QC = ab—#? = —1 (mod 4). 

If f is improperly primitive, % is also improperly primitive since 
it is primitive. 

If ys is properly primitive it represents properly an integer a, prime 
to 2Q0AC. If % is improperly primitive it represents properly an 
integer a, = 2« such that a is prime to 2QAC. Hence ¢ is equivalent 
to a form %, with 4, as the leading coefficient. Transformation 
carrying % into %, carries f into f, = a,2*+.... The contragredient 
(defined by the reciprocal matrix of its transpose) transformation, 


which carries F into the reciprocal F, of f,, leaves C unaltered.* 
We take C, = C. 


6. Characters of f. We may assume, by Lemma 3, that a and 
C or, in case f is improperly primitive, }a and C are prime to 2QA 
and to one another. Multiplying both members of (4.11) by C and 
completing squares, we obtain 
OF = (Cz+ Ry+Sx)?+ (AC—S?)a?+ 2(CT— RS)axy+ (BC— R*)y?. 

(6.1) 
In view of 


AC—S? = Ab, RS—CT = tA, BC— R? = Aa, 


* See, for example, Dickson (4), 16. 
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equation (6.1) may be written as 
CF = X?+ Aba?— 2Atry+ Aay?, X = Cz+ Ry+Szx. (6.3) 

Relations (2.5) and (6.3) may be used to enumerate all the generic 
characters of f. We shall here employ only the latter of the two to 
study characters of f relative to prime divisors of A, and, in case f is 
improperly primitive, also characters with respect to 4 and 8. 

Let M be represented by F. Then (6.3) holds with F replaced by 
M. If p is an odd prime dividing A, then for every M prime to p, 
MC = X? (mod p). Hence 

(M |p) (C\p), 
and f possesses character (F |p). If f is improperly primitive, A is 
even and so are a and 6. Let M be odd. Then, by (6.3), 
CM = X? 1 (mod 4), (6.4) 
and therefore f has the character (—1)!"-». If A is divisible by 4, 
the congruence (6.4) holds modulo 8 and hence f has also the character 
(— 1)8"-)), 
7. Conditions (1.3), (1.5), and (1.6). Let p be an odd prime 
dividing A. Let f= pm. Then 
CX?+0QY? = 0 (mod p) 
in (4.12). If (C|p) = —(—Q|p), X Y = 0 (mod p), and upon divi- 
sion by p (4.12) yields 
Cm = QAz? (mod p), 
whereas f represents all multiples of p. Therefore (C|p) = (—Q!|p). 
But (F'\p) = (C|\p). Hence (1.3) and (1.52) hold. 

[f f is improperly primitive, 

QC = ab—t? = —1 (mod 4), (7.11) 
since Q and C are odd. But (—1)*%-) = (—1)#©-», and hence 
(1.51) holds. Next let A be divisible by 4. In order that f should 
represent all even integers, aCf must represent all multiples of 4 and 
hence all multiples of 8 modulo any power of 2. 

If A is divisible by 8, the only residues of aC modulo 16 are those 
of CX?+QOY? by (4.12). If 

C = —Q-+4 (mod 8), (7.21) 


aCf #8 (mod 16). For if CX?+QY? = 8 (mod 16) and X and Y 
are odd, then 0 = CX?+QY? = C+ Q (mod8), which contradicts 
(7.21). Hence X = Y=0(mod2). Write X = 2X,, Y = 2Y;. 
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Then CX?+QY? = 2 (mod 4), which is impossible since, by 
“ 11), CX?+OY? = C(X?— Y?) (mod4). Hence, in view of 
7.11), C = —Q (mod 8). 

If A = 4 (mod 8), the only residues of QAaz? modulo 64 are 0, 
QAa, 32, where QAa = 8 (mod 16). If (7.21) holds, we show as above 
that CX?+OY? = 8, 24, 32, 40, 56 (mod 64). Then the only possible 
residues modulo 64 which are = 0 (mod 8) are 0, 16, 48, and those 
obtained by adding to them 32 and QAa. Hence the residue 32+-QAa 
is missing. Therefore, as before, C = —Q (mod 8) and (1.6) holds. 


8. Sufficiency of the conditions (1.2), (1.3), and (1.4)—(1.6). 
Since a can be taken as any integer represented by the properly or 
improperly primitive binary form % = ax?+ 2try+by?, we may 
assume that a is (i) an odd prime or (ii) the double of an odd prime 
not dividing QAa according as f is (i) properly or (ii) improperly 
primitive.* Here, as before, C is prime to 2QA. Therefore, by 
Lemma 1, f represents (i) all integers, or (ii) all even integers, if and 
only if O(X,y,2) = X?+QCy2—20 Ryz+QBz* (8.11) 
represents (i) all multiples of a, or (ii) all even multiples of a. Here, 
as before, since C is prime to 2QA, 

(F\p) = (C|p) (8.121) 
for every odd prime p By A and, if f is improperly primitive, also 


(- 1)«F-) pain (—1)e-1 (—1)*-D (—1)KC*-0, (8.122) 


By (8.121) and (1.3) in the case (i), and (1.52) in the case (ii), —QC 
is a quadratic residue of every odd prime divi isor p of A = QO7A = D. 


Case (i). Here f is properly primitive, whence D is not divisible 
by 4, by (1.21), and the above implies that —QC is a quadratic 
residue of D. By (2.3), —QC is also a quadratic residue of a. Hence, 
since a is odd, —QC is a quadratic residue of d = Da and we can 
apply Lemma 2 with L=QC, M = —QR, N =QB. Therefore 
form (8.11) is equivalent to 

®, = 2xy+2*—d2’*. (8.21) 
For z= 0, x = 1 or 2, ®, = 2y+1 or 4(y+1). Next, if d is even, 
whence d = 2 (mod 4), take z = 1, x = 2. Then O, = 4(y+1)+d. 
If d is odd, take z = 1, x = 1, and y = 2u or 2u+1 according as 


* Weber (9). 
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d = 4h—1 or 4h+1. Then ®, = 4(u—h)+2. In both cases ®, is 
equal to any 4/+-2 by choice of y and wu respectively. Thus ®, and 
hence form (8.11) represent all integers. Therefore (8.11) represents 
all multiples of a, and hence f is universal. This completes the proof 
of Theorem 1. 

Case (ii). In this case a = 2q, where q is an odd prime. Write 
D=A = #A,, where A, is odd. Since D is even, k > 1. The deter- 
minant d of the auxiliary form (8.11) is Da = 2*+4A,q. As before, 
by (2.3), —QC is a quadratic residue of a and hence of g. By (1.52), 
—QC is a residue of every odd prime factor p of A,. Hence —QC 
is a quadratic residue of A,qg. It remains to show that —QC is a 
quadratic residue of 2*+1, 

If k = 1, (1.51) and (8.1221) imply that C = —Q (mod 4), whence 
—QC = 1 (mod 4), and therefore —QC is a quadratic residue of 
4 = Qk, 

If k > 2, it suffices to prove that —QC is a quadratic residue of 8. 
As in the preceding case, 

—QC = 1 (mod 4). (8.31) 


Since Q = +1, (1.6) and (8.1222) imply that C = +Q and hence 


—QC = +1 (mod8). Therefore —QC = 1 (mod8) by (8.31) and 
—QC is a quadratic residue of 8. 

Thus —QC is a quadratic residue of d = 2*+1A,q, and we can 
apply Lemma 2 with L=QC, M = —QR, N = QB. . Therefore 
form (8.11) is equivalent to 


2yz+22—dz2?. (8.32) 


Take x = 0 and z = 2. Then (8.32) becomes 4(y+1) which can be 
identified with any multiple 4/ of 4 by choice of y. Therefore (8.32) 
represents all multiples of 4 and hence all even multiples of a. Then 
by Lemma 1, f = 2f, represents all even integers. This proves 


‘Theorem 2. 


9. Universal ternary forms are zero forms. A quadratic form 
is said to be a zero form if and only if it represents zero for the 
values of the variables not all zero. When one compares Theorems 
1 and 2 with the general criterion of Smith* for the existence of 
non-zero solutions of f = 0, one arrives at the following: 


* Smith (8), 410. 
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THEOREM 3. Every universal, classic or non-classic, ternary quad- 
ratic form is a zero form. 

This theorem was first conjectured and proved by Dickson.* We 
shall give yet another, an elementary, proof of this theorem. As 
before, in place of a universal non-classic form f,, we shall consider 
the improperly primitive classic form f = 2f, representing all even 
integers. 

Write A = A, A3, where A} is the largest square dividing A. If f is 
universal, it represents —A,. Since —A, has no square factors, f 
represents it properly. Hence f is equivalent to a form witha = —A,. 
[f f is improperly primitive and represents all even integers, it repre- 
sents properly —A, if this latter is even, and —4A, if —A, is odd. 
For, if —A, is odd, then —2A, is free from square factors and is 
represented, and hence properly represented, by the universal non- 
classic form f, = $f. Therefore, if f is improperly primitive, we may 
take a = —A, or —4A,. Thus in both cases (4.12) becomes 

aCf = CX?+Q(Y?—o%2z?) = G(X, Y,2z), (9.1) 
where « = A, A, or 2A, A,. We see at once that X = 0,y=a,z= 1 
satisfy the equation G = 0. Therefore by (4.13) there exist rational 
and therefore integral solutions not all zero of the equation f = 0. 
To complete the proof it suffices to note that f, is a zero form if 
and only if f = 2f, is such. 

Although the comparison of the Theorems | and 2 and the above- 
mentioned criterion of Smith shows that there are zero forms which 
are not universal, it permits us to state the following theorem for 
the case of determinants free from square factors. 

THEOREM 4. Every zero, ternary, quadratic form f of quadratfrei 
determinant is universal if it is properly primitive and represents all 
even integers if improperly primitive. 

We can also prove it simply as follows. Since f is a zero form it 
is equivalent to a form f, with the leading coefficient zero. Write 

fy, = by? +c2?+ 2ryz+ 2sxz-+- tay, 


gs,, t = gt,, g = (s,t). There exist integers o and 7 such that 


7—8s,o = 1. Transformation y, = t,y+8,2, 2, = oy+7z carries f, 
into a form f’ with a’ = s’ = 0. Hence, omitting the accents, 


f ~ by? +2? + 2ryz+- tay. 


* Dickson (4), 17. 
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The determinant of the latter form is D = —#*c. Since D is free from 
square factors, we may take ¢ = 1.* Thus 
f ~ by?— Dz? + 2ryz+ 2ay. 

Replace x by x+ky—rz. Then, by choice of k, f is equivalent to either 

2ey+y?—Dz* or 2xy—Dz2* (9.2) 
according as b is odd or even, and hence is universal if properly 
primitive. If f is improperly primitive, ) = D = 0 (mod 2) and f is 
equivalent to (9.22), and hence represents all even integers. 


* ifé —1, we change the signs of 2 and 2z. 
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LAPLACE’S OPERATOR IN CURVILINEAR 
COORDINATES 
By T. W. CHAUNDY (Oxford) 

[Received 1 May 1932] 
[In this note I give, by methods of direct differentiation, a simple 
transformation of Laplace’s operator V? to general curvilinear co- 
ordinates u, v, w. Let us denote partial derivatives and partial 
differential operators by suffixes, thus w,, @,,, and let ¥ always denote 
summation over x, y, z. Then we may start with the transformation 
partially effected in the form* 

V2 = (A, B,C, F, G, H) (8,4, &p, Oy)? +(V2u)2,,+(V2v)2,+(V2w)2,,. (1) 

where A=) uv, F =) v,w,, ete. 
To evaluate V?u we regard the three equations 

ut wu + uw =A, 

U,Vz+ Uy Vy+u,v, = H, 

U,W,+U, Wy +U,w, = G, 
as linear equations in w,, u,, u,. So solving them we get 
de A O(v, w) H O(w, u) G o(u, v) | 
eS Oy,z) J Ay,z) | J Ay,z) 
where J = ate il 
O(x, y, Z) 

Differentiate these equations (2) in x, y, z respectively and add to 
give V?u. In this sum the coefficient of A/J will be 

Vey Way) | vy w, | 

> |v, w, | > \O.¢ Wezel 

which vanishes, since the determinants cancel in pairs. Thus A/.J, 
H/J, G/J do not appear in undifferentiated form in the expression 


for V?u, and we are left with 


= O(v,w) 0 [A)\ | O(w, u) oe H = (u,v) @ G 
bi Z Oly, z) Ox (5)+ Z a(y, z) = (3) z O(y, 2) Ox (5). 


O(v,w Ox 
3ut o( hints 2 dl 
Ay,z)  * eu 


* Compare, for instance, Goursat, Cours d’analyse (1910),i. 156 (69), which is 
the corresponding equation for orthogonal coordinates. 
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and so BS A(v,w) @ — J ox bs ee é 


Oy, 2) a 


Ou 0x Cu 


Thus sce F(5)+a(F 2) 
omnes ba d\n 7) OT) * awl \- 


Substitution in (1) and reduction gives 


ees (4, Pe A Pee ee 
= = —d~d.,.+-— = I+C.. C., +-—Cd,..-- = 
a u J ul E vit a u ' v a ul i ot : 


This formula (3) still holds, if we write 
V2 = (a,b,c, f,g, h)(0,, 2,, 8,)*, 
A = (a,b,c, f,9,h)(u,, Uy, Uz)", ete., 


y 
0A OA 0A 

2H = v,— +0, +2, 
Ou, Cu, ou, 


, etc., 


f, g, h are constants. The 


J 


provided that the coefficients a, b, c¢, 
generalization to many dimensions is obvious. 














